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,— | Abstract 

Partition functions of certain classes of "spin glass" models in statistical physics show 
strong connections to combinatorial graph invariants. Also known as homomorphism 
functions they allow for the representation of many such invariants, for example, the 

1 ^ i number of independent sets of a graph or the number nowhere zero fc-flows. 

Contributing to recent developments on the complexity of partition functions |GGJT09l 
ICCL09] we study the complexity of partition functions with complex values. These func- 
tions are usually determined by a square matrix A and it was shown in |GGJT09] that for 
each real-valued symmetric matrix, the corresponding partition function is either polyno- 
Q\ mial time computable or ^P-hard. 

Extending this result, we give a complete description of the complexity of partition 
functions definable by Hermitian matrices. These can also be classified into polynomial 
time computable and #P-hard ones. Although the criterion for polynomial time com- 
putability is not describable in a single line, we give a clear account of it in terms of 
structures associated with Abelian groups. 
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1 Introduction 

We study the complexity of partition functions (or homomorphism functions) on Hermitian^] 
matrices — a class of partition functions which has been studied recently in [LSQSj- Let 
A £ C mxm be a Hermitian matrix and D £ ]g> mxm a diagonal matrix with positive diagonal. 
Let G = (V, E) be a directed graph with possibly multiple edges. The partition function Za,d 
on G is defined by 

Z A ,d(G) = II 4t(u),<j(v) n Ar(u),<T(tO- (1) 

o-iV^— ►[m] nn£-B t>eV 

It is important to note that the product over edges uv £ E counts multiplicities as well. If the 
matrix A is symmetric, the direction of edges in G does not affect the value of the partition 
function. We may then think of G likewise as an undirected graph. In general, however, the 
partition function is well-defined only on digraphs. We refer to D also as a diagonal matrix of 
vertex weights and usually we omit it from the subscript of Za,d if it is the identity matrix. 
Elements of [to] are called spins and mappings a : V — >■ [to] which assign a spin to every 
vertex of G are configurations. Let Ca denote the set of algebraic numbers and let IR& be the 
algebraic reals. 

In this paper we study the complexity of partition functions of the form Za,d f° r a 
Hermitian matrix A £ C^ xm and a diagonal matrix D £ MY lXTn with positive diagonal. Our 



main result (cf. Theorem 1.2, below) is a complete complexity classification of these functions. 



1.1 Examples. 

There are several well-known examples of graph-invariants representable by partition func- 
tions. Let G = (V,E) be a digraph. The number of independent sets of G is given by Zg(G) 
where S is a 2 x 2 matrix which satisfies Sn = and it is one in all other positions. If A is the 
adjacency matrix of a /c-clique without self-loops, then Za(G) is the number of fc-colorings of 
G. 

It will be convenient now and in the following to consider matrices likewise as weighted 
(directed) graphs such that entry Aij is the weight of edge ij and this value being zero denotes 
absence of an edge. Viewed this way, Za(G) is the number of homomorphisms from G to (the 
digraph represented by) A. 

There is an important connection between partition functions and the Tutte polynomial. 
Let G = (V, E) be a graph. For some subset A C E of edges, let c{A) denote the number 
of components of the graph (V,A), then the Tutte polynomial of G, with parameters q,v is 
defined by 

T(G; q ,v)=Y J Q c{A) v lAl - 

ACE 

If the value q is a positive integer then this is also known as the partition function of the 
q-state Potts model. By expanding the above sum over components of G, it can also be shown 
that, for A = J q + v • I q , we have T(G; q,v) = Za(G). Here, J q denotes the all-ones q x q 
matrix and I q the qx q identity matrix. In fact this is also an exact correspondence since (see 
[FLS07]) q £ N is the only way for which T(G; q, v) = Za(G) for a symmetric real- valued A. 

lr That is, Aij — Aji for all i,j, where a denotes complex conjugate. 
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1 INTRODUCTION 



Flows and Hermitian Matrices. Let us turn to an example motivating partition func- 
tions on Hermitian matrices. We start with a simple example (which can be found e.g. in 
[dlHJ93j), given by the matrices 

"-(-11) M 1 -/")- (2 > 

It is not hard to see that, for connected G, the value Zu,Du(G) equals the characteristic 
function of Eulerian graphs. That is, it is 1 if G is Eulerian and otherwise. We will see a 
proof of this in the more general setting of flows which we will introduce now. 

Usually the represent ability of nowhere zero flows by partition functions is derived via 
the Tutte polynomial which has strong connections to both. Deviating from this, a direct 
derivation of this correspondence has been given in [FLS07 . 

Let (5 be a finite Abelian group and S C <3 any subset. For some oriented graph G = (V, E) 
a mapping cj) : E — ^ S and a v E V let 

uv£E vu£E 

Such an assignment <p is called an S-flow, if d&(v) = holds for every vertex v € V. For some 
oriented graph G let fs,<s(G) denote the number of S- flows (w.r.t. (S) in G. Let (5* denote 
the character group of (5. The functions fs.® relate to partition functions in the following 
way. Define two matrices with indices in &*: a matrix A and a diagonal matrix D given by 

D xx = \<8\- X &n&A x , x , = Y,W)-x(g) (4) 

965 

It is not too hard to show that (see [FLS07J for details) for every oriented graph G we have 
Za,d{G) = fs,<s{G)- If we stipulate that S be closed under inversion then A is real- valued 
and symmetric. The orientation of G is then inessential. Particularly, if S = & \ {0} then 
Za,d{G) is the number of nowhere-zero (3 -flows of G. Since this value depends only on the 
order k = \&\ of (5, this counts nowhere-zero k-flows of G. 

If S is not closed under inversion, then generally A is not merely a real symmetric matrix: 
Let, for example <3 = Z3 and S = {2}. Define £ = exp(27ri/3), and consider the character 
group ©* given by x o = (1,1,1), Xi = (1,C,C 2 ) and \2 = (1,C 2 ,C)- The definition of A is 
given by A id = x;(2) • %j(2) which yields 

/ 1 c 2 c \ 
A= c- 2 1 C 1 ■ 

V c- 1 c 1 J 

1.2 Complexity 

A general study of the complexity of partition functions has developed only quite recently, 
although some instances of these functions can be found among the first problems studied in 
the context of counting complexity. It follows, for example, from a result of Simon |Sim77J that 
computing the number of independent sets in a graph is #P-hard. This has also been derived 
independently by Provan and Ball [PB83J. The #P-completeness of counting fc-colorings for 
all k > 3 is shown by Linial in [Lin86j. 



1 . 2 Complexi ty 
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More general complexity results which have a connection to partition functions have been 
obtained by Jerrum and Sinclair [JS93] who show in particular that the partition function of 
the 2-state Potts model is #P-hard to compute. Related to this, Jaeger, Vertigan and Welsh 
|JVW90| study the complexity of computing the Tutte polynomial. 

Results directly related to our work have strong motivational connections to research 
in constraint satisfaction problems (CSP) (see e.g. [FV981 lBu!06l IBK09] ). As a matter of 
fact partition functions and related objects are instances of weighted generalizations of the 
counting constraint satisfaction problem (^CSP). 

Of particular interest are dichotomy results in the spirit of Schaefer's [Sch78| classical 
result on Boolean CSP: Depending on the structure of the constraints allowed in a CSP 
instance, each of these problems is either in P or NP-complete. That such results are not 
self-evident was shown by Ladner |Lad75| : if NP-hard problems are not polynomial time 
computable, then there is an infinite hierarchy of problems of increasing difficulty between 
these extremes. The situation is analogous with counting problems and it is assumed, for 
example, that counting graph isomorphisms is neither in P not #P-complete (cf. |M at79| ). 

Feder and Vardi |FV98j conjectured that such intermediate complexity is not possible for 
CSPs. And in fact, the analogous conjecture for the #CSP problem has been shown to be 
true very recently by Bulatov [B ul08j (see also the new proof by Dyer and Richerby |DR10| ). 
Note that the counting analog of Schaefer's result has been obtained much earlier by Creignou 
and Hermann [CH96J and several weighted extensions of this have been studied recently (see 
[DGJ091 IbDG+OSI ICLX091 IDGJR,09] V 

The line of dichotomy results for partition functions starts with Dyer and Greenhill [DGOO 
who studied partition functions Za on symmetric {0, l}-matrices. This is a counting analog 
of a result by Hell and Nesetf il [HN90] . It turns out that these functions are polynomial time 
computable only if every component of A is either a complete bipartite graph, or a complete 
non-bipartite graph with a loop at every vertex. In all other cases the problem is $R-hard. 

Bulatov and Grohe |BG05j extended this to a dichotomy for all non-negative symmetric 
matrices. To state this result, we need to introduce the notion of blocks of a matrix A. 
Considering A as a graph, blocks arise naturally as follows: each non-bipartite component 
corresponds to one block and each bipartite one corresponds to two blocks. 

Theorem 1.1 (Bulatov &: Grohe [BG05]). Let A G M™ xm be a non-negative symmetric 
matrix. The following holds. 

(1) If every block of A has rank at most 1 then Za is polynomial time computable. 

(2) If A contains a block of rank at least 2 then Za is if-V-hard to compute. 

The fact that this result essentially identifies rank 1 matrices as the polynomial time cases 
renders the result of Bulatov and Grohe widely applicable. Many results on partition functions 
rely on this theorem, as is the case of Bulatov's result [Bul08] . A further example is the work 
of Dyer, Goldberg and Jerrum [DGJ08J on hypergraph partition functions. Dyer, Goldberg 
and Paterson |DGP07| give a dichotomy for partition functions Za where A is an unweighted 
directed acyclic graph. This has been extended to a dichotomy for all directed A with non- 
negative edge- weights by Cai and Chen |CCQ9| . 

There are at least two further extensions of the result in [BG05J to be mentioned. Gold- 
berg, Grohe, Jerrum and Thurley [GGJ T09] gave a dichotomy theorem for partition functions 
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2 PRELIMINARIES 



Za on any symmetric real-valued matrix A, that is, particularly negative entries in A are al- 
lowed. Furthermore, Cai, Chen and Lu |CCL09| extended this dichotomy yet further to all 
symmetric complex valued matrices A. 

The main result of this paper. Here we will extend the work of [GGJT09] into another 
natural direction by studying partition functions on Hermitian matrices. Our main result is a 
dichotomy for these partition functions. As a byproduct, our proof will give a clean account 
of the algebraic structure underlying the proof of [GG JT09J . 

Theorem 1.2. Let A G C™ xm be a Hermitian matrix and D G IR™ xm a diagonal matrix of 
positive vertex weights. Then one of the following is true 

(1) The function Za,d is #P-hard to compute. 

(2) The function Za,d is polynomial time computable. 

Furthermore, there is an algorithm which, given (A,D) decides, whether (1) or (2) occurs. 

Since this work and [CCL09| are both generalizations of the same result [GGJT09J it turns 
out that the proof structure and the polynomial time cases are rather similar from a high-level 
point of view. However, the details are quite different in all three proofs. 

After some preliminaries presented in the next section, we will give a detailed outline of 
the proof of this theorem in Section [3j The details of the proofs will then be given in the rest 
of the paper starting in Section |4j 

2 Preliminaries 

For two counting problems A and B we write A < B if there is a polynomial time reduction 
from A to B. If A < B and B < A holds, we write A = B. 

Computing with Complex Numbers. Recall that Ca denotes the algebraic numbers 
and Ra = Ca H K. For technical reasons, we will always assume that numbers in Ca are 
given in standard representation in some algebraic extension field Q(#). That is, we consider 
numbers in Q(#) as vectors in a d-dimensional Q-vectorspace, where d is the degree of Q(#) 
over Q. It is well-known that for any set of numbers from Ca we can compute a which 
constitutes the corresponding extension field (cf. [Coh93] p. 181). For further details see also 
the treatment of this issue in [DGJ08I IThu09j . 

2.1 Basic Notation 

For natural numbers a < b we let [a, b] := {a, . . . , 6} denote the set of natural numbers from 
a to 6 and we define [b] := [1,6]. 

Graphs. We will be mainly interested in directed graphs G = (V,E) which may have mul- 
tiple edges. We denote edges by uv G E, to explicate that an edge e has multiplicity p we 
sometimes write e p . We denote the out-degree (in-degree, resp.) of a vertex v by deg^w) 
(deg^(f), resp.). The degree of v is dc(v) = deg^(u) + deg^(f) in contrast to the grade 



2.1 Basic Notation 
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9g(v) = deg^(f) — degg(v). A component in a digraph G = (V,E) will be what is usually 
known as weakly connected component. 

Let G be a digraph and p G N. A graph G' is obtained by p-thickening, if each edge e £ E 
is replaced by p many parallel copies of e. Similarly, G" is obtained from G by p-stretching if 
each edge is replaced by a path of length p — that is, a path on p edges. 

Matrices and Vectors. For an to x n matrix A whenever convenient we separate indices 
by a comma Aij to increase readability. Further, A^* (A*j, resp.) denote the ith row (jth 
column, resp.) of A. For an m x n matrix A and sets / C [m], J C [n] the matrix A/j is 
obtained from A by deleting all rows with indices not occurring in / and all columns not in 
J. If I = [m] (J = [n]) we abbreviate this to A* j (Az>, resp.). 
For every p £ Q we define the matrix A^ by 

A (p )= f (A^y , UAijjLO 
*J ' \ , otherwise. 

By tr(.A) we denote the trace of A and 7 m is the m x m identity matrix. Let A be an 
mx n matrix and A' an m' x n' matrix. The tensor product A® A' oi A and A' is the matrix 
defined as follows. Let the indices be given as pairs (i, k) G [m] x [to/] for the row indices and 
(j, I) G [n] x [n'] for the column indices, then 

(A® A') mm = Aij- A kl . 

With the bijections 7r row : [m] x [m'] [m ■ m'] given by 7r row (z,/c) := [i — 1) ■ m! + k and 
7r co i : [n] x [n 1 ] — > [n-n 1 ] given by ir co \(j, I) := (j—l)-n'+l the pairs translate to ordinary indices. 
But most often it will be convenient to explicitly refer to the indices as pairs. In this way we 
will be able to refer to the tiles of A® A' which are the submatrices {A® A')^^^^ = Aij ■ A' . 
The direct sum of A and A' is the (m + m') x (n + n') matrix A © A' defined by 

Aij , if i G [m], j G [n] 

A^- , if z G [m + 1, m + m'] , j 6 [n + 1 , n + n'] 
, otherwise 

Let a, b G C n be vectors. The scalar product of a and b is (a, b) = a\b\ + . . . + a n b n . The 
Hadamard product is the vector ao b = (a±bi, . . . , a n b n ). The conjugate Hadamard product is 
a • b = a o b. 

We denote by U the set of all complex numbers a of absolute value 1. The algebraic 
numbers of absolute values 1 are denoted by Ua- Similarly, U w denotes all w-th roots of 
unity. A matrix A G \] mxn is called normalized if its first row and column are constantly 
1. The expression A denotes the complex conjugate of A. A matrix A G U mxn is called a 

— T 

complex Hadamard matrix if AA = nl n , that is all of its rows are pairwise orthogonal. Note 
that this implies that the columns of A are pairwise orthogonal as well, which follows from 
elementary group theory. 

Matrices, blocks and their connection to graphs. An mxn matrix A is decomposable, 
if there are non-empty index sets / C [to], J C [n] with (I, J) / [to] x [n] such that A^ = 
for all (i, j) G I x J and all (i, j) G / x J. A matrix is indecomposable if it is not decomposable 
and a block of A is a maximal indecomposable submatrix. 



( A © A% = < 
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2 PRELIMINARIES 



Let Abe an mxm matrix. The graph underlying A is the digraph defined by Ga = ( [m] , E) 
with E := {ij \ Aij ^ 0}. For Hermitian matrices the digraph Ga is symmetric in the sense 
that ij G E 1 iff ji £ E. For considerations on the connectedness of this graph it is therefore 
sufficient to consider it as undirected. 

We call a matrix A connected if Ga is connected. We call a submatrix An a component 
of A, if Cr^pT] is a component of Ga- If GUpT] is non-bipartite, then An is a block of A. Note 
in particular that a connected non-bipartite matrix is is a block, and it is thus identical with 
its underlying block. If otherwise Ga[I\ is bipartite then An has an underlying block B such 
that 




2.2 Generalized Partition Functions 

Let A 6 C mxm be a matrix and D 6 C mxm a diagonal matrix. Let G = (V, E) be some given 
digraph. A pinning of (vertices of) G with respect to A is a mapping (ft : VF — >■ [m] for some 
subset If C ^ of vertices of G. For a given graph G = (V, E) and a pinning 0, we define the 
partition function 

Za,d{(I>,G)= ^ Y\_ A <j{u),<j(v) \\ D a(v),<r(v)- 
^C(7:^[m]««eE ueV\def(0) 

Note in particular, that the sum is over all configurations a : V — > [m] which extend 
the fixed given pinning (ft and, for technical reasons, the terms D C ^ V ^ (T ^ for v G def(0) 
are excluded from the above expression. The weight of the configuration a is the term 
UuveE A a(u),a(v) Uv€V\dci(<p) D a{v),a{v)- Whenever 4> is trivial, that is, def(» = and D is the 
identity we omit the terms D (</>, respectively) in the expression. Let EVAL pm (^4, D) denote 
the problem of computing Za,d{<P^G) on input G,4>. Analogously, EVAL(v4, D) denotes the 
problem restricted to empty input pinnings. 



Congruential Partition Functions. We say that a matrix A £ C™ xm is to-algebraic, if 
for every non-zero entry Aij there is some w-th root of unity ( such that Aij = \Aij\ ■ (. 
Working with partition functions on w-algebraic matrices, it will be convenient to transition 
to a different formulation of such problems. 

Let, for a d £ Z the expression [d] denote the modulus of d after division by oj. As oj is 
defined relative to A, the exact meaning of [d] will always be clear from context. A family 
)cez^ of diagonal mxm matrices will be called a family of congruential vertex 
weights. The congruential partition function on A, and D = (D^) c <zi^ , for a digraph G and 
a pinning (ft is defined by 

zam<p,g)= £ n ^h^) n 

where V' = V \ def(<ft). Recall that d(v) = deg + (t>) — deg~(w) denotes the grade of v in G. 
The connection to ordinary partition functions easily follows from the definitions: 
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Lemma 2.1. Let A G C™ xm be a Hermitian u-algebraic matrix and D G M™ xm a diagonal 
matrix of positive vertex weights. Let D = (Z)I c H) cg ^ such that £H C J = D /or a// c G 

Za,d{<I>, G) = ZA,T>{(j)i G) for all digraphs G and all pinnings <j). 

Correspondingly, EVAL pm (yl, D) denotes the evaluation problem associated with this parti- 
tion function on input (cf>, G). 



3 The Proof of the Main Theorem — An Itinerary 



We will now give a detailed overview of the most important steps in the proof of Theorem 1.2 
Given an arbitrary Hermitian matrix A and diagonal matrix D it is rather routine to transition 
to the case of connected A and the more general problems EVAL pm (^4, D). 

Starting from this point, we will filter the abundance of problems EVAL pin (^4, D) until 
we are left with only the polynomial time computable ones. This process of filtering will 
be performed by transforming in several steps the problem EVAL pm (A, D) into a problem 
satisfying certain conditions. In each step more and more such conditions will be accumulated 
and if any of these conditions fails for a certain problem it will be shown that this is due to 
the problem being #P-hard. 

The first step is to gather conditions such that the resulting problem can be transformed 
into a problem of the form EVAL pm (^4', D) where the block underlying the connected matrix 



A' is a complex Hadamard matrix (cf. Lemma 3.4) — we call such matrices Hadamard 



components. Once we have achieved this, we will be able to gather more conditions on these 



problems (cf. Theorem 3.5). And eventually, we will have to solve the non-trivial task of 
providing a polynomial time algorithm for the remaining problems — Theorem |3.8| yields 
this. 

Some words are in order for the additional statement of Theorem 11.21 It is claimed that 
there is an algorithm deciding whether given matrices A, D give rise to polynomial time 
computable partition functions or #P-hard ones. This not very surprising fact will follow 
straightforwardly from the proof of Theorem \1.2\ We give details on the decision algorithm 



only if they do not follow straightforwardly from the details of the proofs given. 



3.1 Prom the General Case to Connected Matrices 

We start by transforming the original problem EVAL(j4,D) into an equivalent problem 
EVAL pm (yl / , D) with A' an w-algebraic matrix. The transition to w-algebraic matrices is 



captured by the following lemma. Its proof can be found in Section 4.3 



Lemma 3.1 (The Arithmetical Structure Lemma). Let A G C™ xm be a Hermitian 
matrix and D E M™ xm a diagonal matrix of positive vertex weights. There is an u G N 
and an co-algebraic matrix A' whose entries satisfy l-A^I G N and 

EVAL(A',D) = EVAL(A, D). 



With the following lemma we then take the step which allows us to pin certain vertices of the 
input. The proof will be given in Section 4.4 
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3 THE PROOF OF THE MAIN THEOREM — AN ITINERARY 



Lemma 3.2 (Pinning Lemma). Let A G C™ xm be a Hermitian matrix and D G M™ xm a 
diagonal matrix of positive vertex weights. Then 

EVAL p/n (A, D) = EVAL(A D). 



The transition to congruential partition functions is a consequence of Lemma 2.1 
Corollary 3.3. We have EVAL P '"(A, D) = EVAL pin (A, D). 

3.2 Prom Connected Matrices to Hadamard Components 

We start under the following preconditions. We are given an w-algebraic Hermitian matrix 
A G (T^mxrn w j 1 i c ] 1 satisfies that \A{j\ G N for all i,j G [m] however this being in N is not 
essential. We say that a matrix H and a family of diagonal matrices D define an (H— STD) 
problem, if the following holds 

(H STD) There are n, oj G N and D = (£H c I) cg z w is a family of diagonal matrices such 
that 

• H G U" xn is a normalized Hermitian Hadamard matrix.. 



• We have L>M = j n> £)H G ({o} u U A ) nxn and = flW for all c G Z w . 

Further, an evaluation problem EVAL pm (74, D) is called a (B— H— STD)-problem if it satisfies 
the following conditions. 

(B H STD) There are n,u G N such that A is a connected Hermitian matrix with under- 
lying block H and D = (L>M) c6Zw is a family of diagonal matrices. Further, 

• H G U" xn is a normalized complex Hadamard matrix. 

• For all c G 1 W there are diagonal matrices D^' TOW , £>H> co1 g ({0} UU A ) nxn such 
that 

£)14 _ £)[c],row £j[c],col_ 



• DM = / 2n , and Z?I- C 1 = DH for all c G Z u . 

The following lemma shows how we can reduce the complexity classification of evaluation 
problems on arbitrary connected Hermitian matrices to those on Hadamard components. Its 
proof will be given in Section [5] 

Lemma 3.4. Let A be a connected co-algebraic Hermitian matrix and D a diagonal matrix 
of positive vertex weights. Then either EYAL pm (A, D) is ^V-hard or the following holds. 

(1) If A is not bipartite then there is an u-algebraic Hermitian matrix H and a family of 
diagonal matrices D = (D^)^^ which define an (H STD) problem such that 

EVAL p/ "(^, D) = EVAL p '"(i7, D). 

(2) If A is bipartite then there is an to-algebraic Hermitian matrix A 1 and a family of 
diagonal matrices D = (D^)^^ which define a (B H STD) problem such that 

EVAL pin (A,D) = EVAL^^^D). 



3.3 Hadamard Components 
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3.3 Hadamard Components 



After applying Lemma 3.4 we are left with the task of classifying the complexity of (H— STD) 
and (B H STD) problems. We will now develop further conditions on these problems which 
will eventually lead us to polynomial time computable problems. The definition of these 
conditions will be given in the following. We start with the statement of the theorem whose 
proof can be found in Section [6j 

Theorem 3.5. Let A be an uj -algebraic matrix and D a family of diagonal matrices defining 
either an (H-STD) problem or a (B H STD) problem. Then either EVAL p/ "(A,D) is 
jfV-hard or the following holds. 

(1) If A, D define an (H STD) problem then they satisfy conditions (GC), (Rl) through 
(R5) and the Affinity Condition (AF). 

(2) If A, D define a (B H STD) problem then they satisfy conditions (GC), (B Rl) 
through (B R5) and the Affinity Condition (B— AF). 

The first condition (GC) stipulates that the Hadamard matrix H underlying every (H— STD) 
and (B H STD) problem has a certain group structure. 



The Group Structure of H. Let R{H) := {Hi^ \ i G [n]} denote the set of rows of H. 
The Group Condition is defined as follows. 

(GC) H it * o H jt , G R(H) and fl£ o fl£ G R(H T ) for all i,j G [n\. 

The name "Group Condition" is justified as it guarantees that R{H) and R(H T ) are Abelian 
groups under the Hadamard product o. To see this, note that the definition of (GC) is known 
as the subgroup criterion for finite groups. 

To introduce all further conditions, we will consider the non-bipartite case of (H— STD) 
problems and the bipartite (B H STD) problems separately. Let us begin with the non- 
bipartite case. 



A Group Theoretic Account of (H STD)-problems. When we are working with 
Abelian groups (5, unless defined otherwise, the group operation will be denoted by + and 
the neutral element will be 0. Let an w-algebraic matrix H and a family D = (-DW)z^ define 
an (H STD) problem. Further, assume that H satisfies the group condition (GC). 

Fix the canonical isomorphism & u : 7L U — > U w where U w denotes the set of w-th roots of 
unity. That is, we have X^{a) = exp(^^). By the group condition (GC) the set R(H) of 
rows of H is an Abelian group under the Hadamard product. Further by the Fundamental 
Theorem of Finitely Generated Abelian Groups R(H) is isomorphic to some Abelian group 
(5 = Z qi © • • • © Zq z with qi some prime powers. 

Note first that the order of each element g G divides u. To see this, recall that g G (3 
corresponds to a row Hi^ by the above mentioned isomorphism. The order of Ui^ is the least 
common multiple of the orders of its entries Hij. However, each of these entries is an cj-th 
root of unity. 

Let / : R(H) — > (3 be an isomorphism. We see particularly that for all i G [n] there is a 
unique a £ & such that /(iJ^*) = a and — as H is Hermitian — /(i?*^) = f(Hi^) = —a. We 
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may therefore relabel the entries of H by elements of such that H a ^ o H^ = H a+ b,* for all 
a,b £ 0. In this way, if for i,j £ [n] we let a, b £ be their corresponding group elements, 
then we have Hij = -ff a ,-f>- We thus obtain H a _b = Hij = Hji = H^-a- 
We further define an operator (*, *) w : x — > Z w such that 

H a _ 6 = jr w [(a, 6) w ] for all a,b £ 0. 

It thus follows directly from our considerations that 

Lemma 3.6. XTie operator (*, *) w /ios i/ie following properties. 

• For all g £ & the mappings (g, *) w and group homomorphisms. 

• For a// g, h £ <5 we have (g, h) w = — (/i, g)^. 

We summarize these properties by saying that the operator (*, *) w is skew-bilinear. Extending 
this description, we would like to describe the non-zero entries of the diagonal matrices 
in a similar way, although our reasoning so far does not allow us to stipulate that Z)j£f / 
implies D{ c } a £ U w . 

Therefore we extend Z w to some set $7 and SE W to some mapping X : fl — > U such that 
for all c £ Z u and all a G 0, if dJ^J / then there is some a £ Q such that 3C\a\ = D^f a . 
We may further assume that £1 is itself an Abelian group and is a homomorphism. (Note 
that this is possible if we just define as the real interval [0,lj) and JT(a) = exp(^ • a)). 
Altogether, from conditions (H— STD) and (GC) we have derived the following. 

(Rl) There is an Abelian group f2 D Z w and a homomorphism : 17 — > U such that |z„ 
is the canonical isomorphism between the additive group Z w and the multiplicative 
group U w . 

(R2) R(H) is isomorphic to some Abelian group = Z 9l © ... © Z 9z for q±, . . . , q z some 
prime powers. There is a skew-bilinear operator (*, *) w : x (5 — ► Z w such that 

H a _ b = 3£ [(a, b) u ] for all a, b £ 0. 

Furthermore, the order of each g £ <8 divides w. 

We shall define three further conditions. The first of which will be shown to hold in the 
following. 

(R3) For every c £ Z w let A c = {a £ | dI% / 0} be the support of F>M. Then there is 
some f3 c £ and a subgroup C of such that 

Ac = P c + C . 

(R4) For every c£Z u with Z)M / there is a mapping p c : C — > $7 such that 

= ^ [Pc(a - &)] for all a £ A c . 



(R5) For every c£Z u with I)H / Owe have p c (0) = 0. 



3.3 Hadamard Components 
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Note that, by the above, condition (R4) holds once we made sure that (R3) is satisfied. We 
will define a further criterion on the functions p c describing the non-zero vertex weights. We 
call this criterion the Affinity Condition for D' c ': 

(AF) For all cG Z u there is a mapping 7 C : © c — > © such that following is true. We have 
p c (y + x) - Pc(x) - p c (y) = ('Yc{y),x)«, for all x,y G © c . (5) 



A Note on Decidability. Apart from the conditions presented in other sections, the 
decidability of the conditions given here may, at this point, be a bit obscure to the reader. 
We shall give a bit of insight into this, now. Assume, we are given a matrix H and a family 
D of diagonal matrices defining an (H— STD) proble m an d we want to decide whether this 

it turns out that EVAL pin (if, D) 



3.5 



problem is polynomial time computable. By Theorem 
is #P-hard unless all of the conditions just described are satisfied. To be precise, as we have 
not given an algorithm which can compute the representation as given in (Rl) to (R5), for 
the cases of #P-hardness Theorem |3.5| does only prove the existence of such a reduction. 
Eventually however, if the evaluation problem is polynomial time computable, it will turn 



out that each of the matrices D' c ' satisfies the following 

For all c£Z u all non-zero entries of are in U2 W . (6) 

In a first step we may thus check whether this is the case. If so it follows straightforwardly 
that the Abelian group defined in condition (Rl) is finite and all of the conditions defined 
can be computed straightforwardly. 



A Group Theoretic Account of (B H STD)-problems. Let A be an w-algebraic 
matrix and D a family of diagonal matrices define a (B— H— STD)-problem. Assume further 
that the block H underlying A satisfies the group condition (GC). 

Fix the isomorphism : Z w — > as before. By the group condition (GC), the set 
R{H) of rows of H and the set R(H T ) of its columns are Abelian groups under the Hadamard 
product. Hence, by the Fundamental Theorem of Finitely Generated Abelian Groups R{H) 
is isomorphic to some Abelian group © r0 w = Z 9l @ ■ ■ ■ @ 1^ qz with qi some prime powers. 
Further, R(H T ) is isomorphic to some Abelian group © co i = Z ? / © • • • © Z q / ; with q\ some 
prime powers. 

Note first that the order of each element g G © r0 w (<?' G © co i resp.) divides u. To see this 
recall that g G (5 r ow corresponds to a row Hi t * by the above mentioned isomorphism. The 
order of Hi^ is the least common multiple of the orders of its entries Hij and each of these 
entries is a w-th root of unity. 

Let / row : R(H) — > © row and f co \ : R(H T ) — > <3 co i be group isomorphisms. We see 
particularly that for all i G [n] there is a unique a G <3 r ow such that / r0 w(^,*) = o, and the 
analogue of this holds for the columns. We may therefore relabel the entries of H by elements 
of ©row and ©coi such that H a ^ o iJ 6 * = H a+ b^ for all a, b G © row and ff* jC o H*^ = -ff*, c +d 
for all c, d G © co i- We further define an operator (*, *) w : © row x ©col — > Z w such that 

H a ,b = ^lj [{a, b) u ] for all a G © row , b G © co i. 



2 The interested reader may find the proof of this in Lemma 



7.4 
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Be aware that, although we denote this operator in the same way as the corresponding 
one in the non-bipartite case, the two operators are different. However, it will always be 
unambiguously clear from the context, which of the two operators we are referring to. 

We say that the operator (*, *) w is bilinear if the following holds: For all g G <3 row and 
all g' G ©col the mappings (g, *) w and (*,g')uj are group homomorphisms. It follows directly 
from the above considerations that 

Lemma 3.7. The operator (*, *) w is bilinear. 

We would like to describe the non-zero entries of the diagonal matrices Z)M in a similar way, 
although our reasoning so far does not allow us to stipulate that D^ a / implies G U w . 
Therefore we extend Z^ to some set f2 and to some mapping S£ : f2 — > U such that for 
all c£Z u and all a G [2n], if D^l / then there is some a G such that JT[o;] = We 
further assume that O is itself an Abelian group and is a homomorphism. 

Altogether, from conditions (B-H-STD) and (GC) we have derived the following. 

(B Rl) There is an Abelian group and a homomorphism : 17 — > U such that \z u 
is the canonical isomorphism between the additive group Z w and the multiplicative 
group ILL. 

(B-R2) R(H) is isomorphic to some Abelian group © row = Z 9l ®. . .©Z 9z for gi, . . . , q z some 
prime powers. R(H T ) is isomorphic to some Abelian group (5 co i = Z ? / © . . . ©Z ? / / for 
<7i, ■ ■ ■ ,q' z i some prime powers. There is a bilinear operator (*, *) w : (5 row x ©col Z w 
such that 

= X [(a, 6) w ] for all a G © row , b G <3 C ol- 
Furthermore, the order of each g G <3 row and each g' G <5 C oi divides uj. 

We shall define three further conditions. The first of which will be shown to hold in the 
following. 

(B-R3) For every c G Z w let A c , row = {a G © row | d[%' row / 0} be the support of D'*™ 
and let A CjCol = {a G co1 | £>f|' co1 / 0} be the support of D^' co \ There are 
/3 c ,row G © r ow, Pc,co\ G ©col and subgroups <5 Cir0 w of © row and CjCO i of (5 co1 such that 

Ac,row = /?c,row ~i~ ©c,row and A C C ol = Pc,co\ ©c,col- 

(B R4) For every c G Z w there are mappings p c ,row '■ ©c,row — > ^ and p c , co \ '■ ©c,col — > & 
such that 

D l4,ro W = % [pc mw{a _ /3c row) ] for all a G Acrow . 
^W.coi = ^ [pccoi(a _ for all fl G Accoi 

(B-R5) For every c G Z w the following holds. If Z)W' row / then p c , ro „(0) = and if 
L>H' co1 / then /9 c , co i(0) = 0. 

Note that, by the above, condition (B R4) holds once we made sure that (B R3) is satisfied. 
As in the description of the (H-STD) problems we will give a further condition on the 
mappings p c ,row, Pc.col describing the non-zero vertex weights. Again, we call this the Affinity 
Condition for Z)M, which reads 



3.4 The Polynomial Time Case 
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(B-AF) For all cGZ u there are mappings 7 Cirow : & c<mw -> © co i and 7 CjCO i : <5 CjCO i -> <5 row 
such that the following is true. 

Pc,mw 

(y + x) — Pc,row\%) Pc,row (y) = ( •E) Tc,row for all x, y £ ® c ,row (7) 

Pc,co\(y + x) - p c ,coi(x) - p c ,co\{y) = (7c,col(y) ! aj>w for all sc, y G <3 c , C oi- (8) 

Decidability. Note that the question whether a given (B H STD) problem is polynomial 
time computable or #P-hard can be solved analogously to the non-bipartite case. Recall the 
discussion given there. 

3.4 The Polynomial Time Case 

By now we have arrived at (H— STD) problems and (B H STD) problems satisfying all 
conditions we have defined so far. It remains to show that this implies polynomial time 
computability of the corresponding problems. 

Theorem 3.8. Let A be an oj-algebraic matrix and D a family of diagonal matrices. The 
problem EVAL p '"(j4, D) is polynomial time computable if the following holds. 

(1) If A, D define an (H STD) problem then they satisfy conditions (GC), (Rl) through 
(R5) and the Affinity Condition (AF). 

(2) If A, D define a (B H STD) problem then they satisfy conditions (GC), (B Rl) 
through (B R5) and the Affinity Condition (B— AF). 

This theorem will follow by reducing the partition functions satisfying the above conditions 
to a problem which has been shown to be polynomial time computable in [CCL09J. The proof 
details will be given in Section [7| 

3.5 The Proof of Theorem Q 

With the high-level results we have derived in the previous sections we are now able to prove 
the Main Theorem |1.2| In a first step, let us show how we can easily derive the result for 
problems EVAL pm (A, D) on connected matrices A. 

Lemma 3.9. Let A £ C™ xm be a connected u-algebraic Hermitian matrix and D G M™ xm 
a diagonal matrix of positive vertex weights. Then either EVAL p '"(t4, D) is ftP-hard or poly- 
nomial time computable. 

Proof. Let us consider the case first that A is not bipartite - the bipartite case follows by 



analogous reasoning. Assume that EVAL pm (^4, D) is not #P-hard. Then Lemma 3.4 implies 
that there is an w-algebraic Hermitian matrix H and a family D = (Z)W) c6 2 w of diagonal 
matrices which define an (H— STD) problem such that 

EVALP in (^l,D) = EVALP in (#,D). (9) 



Therefore EVAL pm (ff, D) cannot be #P-hard. By Theorem 3.5 the matrix H and the family 



D thus satisf y co nditions (GC),(R1) through (R5) and the Affinity Condition (AF). Hence 



by Theorem 3.8 the problem EVAL pin (i7, D) is polynomial time computable. Therefore, 



equation Q implies that EVAL pin (A, D) is also polynomial time computable. n 
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Proof (of Theorem\Tjfy. Let A G C™ xm be a Hermitian matrix and D £ M™ xm a diagonal 
matrix of positive vertex weights. Assume first that A is not w-algebraic for some u) E N, then 
the Arithmetical Structure Lemma |3 . 1 1 implies that there is some w 6 N and an w-algebraic 
matrix A' such that EVAL(^4, D) = ~EVAL(A' , D). We therefore assume for simplicity that A 
is u- algebraic. We have 

EVAL pin (A D) = EVAL(A, D) (10) 



by the Pinning Lemma 3.2 Let I be a partition of [m] into subsets I £ X such that each An 



is a component of A. The following is true for each / 6 X. 

C/aim 1. If EVAL pin (y4//, D//) is #P-hard then EVAL(^,L>) is #P-hard. 

To see this, let G, <f> be the input to EVAL p,n (Au, Djj) for some J £ I. We may assume 
that G is connected and def(^) / 0. Therefore, we have Z Aii ^d i1 {4>^G) = Z A £,(4>,G) and 
thus ENALV m (A n ,D H ) < EVAL pin ( J 4, D). By equation ([To]) this proves Claim[l) We further 
have, for all digraphs G, 

Z a ,d(G) = ^ Z Aii ,d u {G). 
/ex 

This proves that EVAL(j4,D) is polynomial time computable if for all I £ X the problem 



EVAL(Ajj, Du) is. By Lemma 3.9 we see that for each component Ajj of A either the 
preconditions of Claim [l] are given or EVAL pm (j4//, Djj) is polynomial time computable. 
This finishes the proof. n 



4 Basic Technical Results 



In this section we will give proofs of many technical results we will need later on. Most 
importantly, we will prove the Arithmetical Structure Lemma 3.1 in Section 4.3 and the 
Pinning Lemma 3.2 in Section 4.4. The proofs of these results require a considerable amount 
of other tools which we will develop in the following. 



4.1 General Principles 

Most reductions which we will present here have to deal with arbitrary input digraphs G and 
pinnings (j). To avoid tedious reasoning on trivial or vacuous cases we will introduce some 
general principles which we will apply whenever necessary. 

Lemma 4.1 (Connectedness Principle). Fix one of EVAL P '"(A, D), EVAL p/n (A D) or 
EVAL(yl, D). If we know how to solve the problem in polynomial time given that the input is 
restricted to connected digraphs, then the general problem can be solved in polynomial time. 

Proof. We give the proof for EVAL pm (A, D), the other cases are analogous. Let G = (V,E) 
be a given digraph with components G±, . . . , G c . Each pinning 4> of vertices of G to entries of 
A can be partitioned into pinnings (possibly empty) cpi, . . . , <p c corresponding to the vertices 
of the components of G. Then 

c 

Z A>D ( ( j ) ,G) = Y[Z^D^,Gi). 

i=i 



Since we know how to compute the values Z At n(<t>i, Gi) in polynomial time, the value Z Aj d((J), G) 
can be computed by the above equation. n 



4.2 Technical Lemmas 
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Lemma 4.2 (The Pinning Principle). Fix one of EYAL pin (A, D) or EVAL P '"(A, D). As- 
sume that this problem can be solved in polynomial time if the input allows only non-trivial 
pinnings. Then the general problem can be solved in polynomial time. 

Proof. We give the proof for EVAL pm (A, D), the other cases are analogous. Assume that A 
is an m x m matrix and let G, <f> be an instance of EVAL pm (A, D). If def(0) = 0, fix an 
arbitrary vertex v £ V and let (pi for every i £ [m] be the mapping defined by v t— >• i. We 
have 

m 

i=i 

Each value Z^£>((f>i, G) can be computed in polynomial time. n 

For A an m x m matrix and ir : [m] —> [m] a permutation, define A wn by (A nn )ij = 
A^u\^u\ for all i,j £ [m]. It is straightforward to see that simultaneous permutation of the 
rows and columns of A do not alter the partition function: 

Lemma 4.3 (Permutability Principle). Let A,D 6 C™ xm and tt : [m] — > [m] a permu- 
tation. Then EVAL p '"(t4, D) = EVAL p '"(A 7r7r , D nn ). This holds analogously for congruential 
partition functions. 



4.2 Technical Lemmas 

Lemma 4.4. Let < x\ < . . . < x n be positive in ci, . . . , c n € Ca- There is a po G N 
such that for all p > po, the following equation holds if, and only if, all coefficients Ci are 
zero. 

n 

o = 5><- ( u ) 

i=l 

Further po is computable from input xi, . . . , x n , c\, . . . , c n . 



Proof. Note that, if = c\ = . . . = c n , then equation (11) is satisfied for all p > 0. Assume 
therefore that there is at least one non-zero q. We perform induction on the maximum index 



m < n such that c m ^ 0. If m 
For m > 1 note that 



1 then Y^i=i c i x i- = °i x i an d the proof follows for po = 0. 



E 

i=l 



C j .1 ■ 



> 



Cm.Xf 



m—1 

E 

i=l 



( j J ■ 



and the right hand side is larger than zero, if |c m |xm > YULi^ \ c i\ xP i- If an c±, . . . , c m _i are 
zero, we are done with po = and let x = maxjxi | c« 7^ 0, i £ [m — 1]} otherwise. We see 



that equation (11) is not satisfied if 



x m \P 



> 



m—1 

E 

i=l 



\Cj\ 
I C-m I 



As x m > x this inequality is satisfied for large enough p. We define po such that this holds. n 
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4.2.1 Twin Reduction. 

Let A e C mxm be Hermitian. We say that two rows and Aj t * are twins, if A^* = Aj^. 
A matrix A is twin-free if it has no twins. Twins induce an equivalence relation on the rows 
of A. Let Ii,...,Ik be the equivalence classes of this relation. Since A is Hermitian these 
classes are also the equivalence classes of the twin relation on the columns of A. The twin 
resolvent of A is the k x k matrix [A], defined by 

[A]ij = A^ tU for some //G/j and v G Ij. 

We say that [A] is obtained from A by twin reduction. The twin resolution mapping r : [m] — > 
[k] of A is defined such that fx £ ^t(h) f° r au I 1 £ [ m \- Hence [A] T ^^ T ^ = Aij for all i,j £ [m\. 

To use twin reductions in the context of partition functions, we need to consider their 
effect on the diagonal matrices D associated to A. We will see that the following diagonal 
k x k matrix captures this effect 

^K ] = E^ foralHG [ fc ]- ( 12 ) 

Analogously, for a family D = (-D" c ") c g^ of congruential vertex weights, we define a family 
Dt A l = (-D^'^)c6Z w °f diagonal k x k matrices by 



Df} M = Of} for all i G [k] and all c e Z u 



Lemma 4.5 (Twin Reduction Lemma). Let A 6 c mXTra 6e Hermitian. Let I\, . . . ,1^ be 

the equivalence classes of the twin relation on A and r be the twin resolution mapping of A. 
The following is true for all digraphs G and pinnings (ft: 

(1) Let D = (£>M 

)ceiu> a family of congruential vertex weights, then 
Za,t>(4>, G) = Z [A]iG[A ] (r o(f>,G). 

(2) Let D be a diagonal m x m matrix of vertex weights. Then 

Za,d{4>, G) = Z [a]>d1 a] (r ocf), G). 

Proof. Let us start by proving (1). For a digraph G = (V, E) and a pinning (ft, let V = 
V\def((ft). We have 

z^G) = £ n Mo,™ n D f$U) 

= E II [A] T o<r(u),TO<r(y) ]^[ D i(i)^( v ) 



where the second equality follows from the definition of r. As for all configurations a : V — > 
[m] we have r o a : V — > [k] , we can partition the a into classes according to their images 



4.2 Technical Lemmas 



19 



under concatenation with r and obtain: 



e e n n 



TOa=a / 



e n [^]^(«)^'(«) 

Totj>Cu':V^-[k] uv&E 



(v) 



. <£C>:F->-[m] «£V' 

\ TO(J=(t' 



J 



Fix some a' : V — > [k]. For a : V — > [m] we have r o a = o 7 if and only if <r' ({*}) = <7 _1 (Xj) 
for all i£ [k]. Define for all i G [fc] the set := cr /_1 ({?}) and the mapping <pi := fdcf(»nVi) 
then 



\^ TT n [S(«)] 
Z> 11 a(u) j( r(<;) 



</>Ccr:V-S-[m] vSV 



' 4 11 cr(y),cr(v) 

tj>C(r:V-t[m] veV 
Vie [A;]: o-CV^C/i 

n r TT z? I0(,,)I 

11 1^1 11 ^ffi(«),<Ti(«) 

i=l faCai-.Vi^U «6V^\def(0i) 

n n e 

i=l ueVi\def(^) i/SA 



n n [A],[e(«)] 



^ev 



Hence 



For the proof of (2) notice that flW = D implies D^I-M = for every c G Z w . The claim 
hence follows by Lemma |2.1| n 

Corollary 4.6. We have EVAL p/n (A, D) = EVAL p/ "(L4], D^). 
4.2.2 Root Of Unity Transfer 



Lemma 4.7 explains a way to transfer certain values from the Hermitian matrix A G iry^xm 
into the vertex weights of the partition function. 

Lemma 4.7 (Root Of Unity Transfer Lemma). Let A G C mxm be Hermitian and u- 
algebraic. Let D = (-D^)ceZa, be a family of congruential vertex weights and let IT be a 
diagonal mx m matrix whose diagonal entries are oj-th roots of unity. Define A 1 = UAH and 
a family of diagonal matrices D' = (D' M ) c&hj by (-D')M = WD^. 

There is a polynomial time computable function /n such that for every digraph G and 
pinning <p we have. 

ZaM^G) = M<j>)-Z A ,p>{<j>,G). 
Furthermore, /n is non-zero for all (p. 
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Proof. Let G = (V, E) be a given digraph and <f> a pinning. With V' = V \ def(0) we have 



<r:V->[m] troS-E v£V 
aD<j> 

r(u)<T("t>) 

cr:V-s>[m] troe-E »6l/' 



TT rr 9 ^ IV TT a' TT n 5W d 19(v)} 



v i;edef(0) / a:V->[m] uvEE veV 

oOfl 



Now IT 9 ^) = nl 9 ^)l as all diagonal entries of II are w-th roots of unity. Therefore, 

\v£def(<f>) / <r:V->[m] uvGE v&V 

n n ^W) ^2 n ^w^) n ^s^V) 

K v£def(4>) J a:V^[m]uv£E v£V' 

o-D<t> 

^Sdef(</>) / 

Note that the factor ^n^edef(0) ^(uj <X*>)) * s nonzero an< ^ polynomial time computable. In 
particular, this factor equals 1 if def(0) = 0. n 

The above lemma directly yields, 

Corollary 4.8. We have EVAL p/ "(A,D) = EVAL P '"(A / , D'). 
4.3 The Arithmetical Structure Lemma 



In this section we will give the proof of Lemma 3.1 We need some preparation 



EVAL(A, D) and COUNT(A, D). Let A G C™ xm be a matrix and G = (V, E) a digraph. 
We define a set of potential weights 

W A {G) := I Yl A 7j* j I Yl m ij = \ E l and > 0, for all i,j G [m] > . (13) 

^*>j'e[m] i,j6[m] J 

For some u; G Ca denote the set of configurations with weight w by 

Ca(G, w):=<a:V->[m]\w= J| A a{u)a[v) \ . 

I uvGE ) 



4.3 The Arithmetical Structure Lemma 
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Define the value 

N A ,d{G,w) := ^ II D cr(v)*(v) 

aeC A (G,w)vev 

COUNT(A, D) denotes the problem of computing Na,d{G,w) for some given digraph G = 
(V, E) and a weight w G Ca- We use an interpolation technique from [DGOOj : 

Lemma 4.9. Let, for some fixed 6 G Ca let x\, . . . ,x n G Q(9) be pairwise different and 
non-zero. Let b\, . . . , b n G Q(9) be arbitrary such that 

n 

bj = °i x l f or a ll 3 ^ H- 
1=1 

Then the coefficients c\, . . . ,c n are uniquely determined and can be computed in polynomial 
time. 

Lemma 4.10. For every matrix A £ C™ xm and D £ C™ xm we have 

EVAL(A, D) = COUNT(A, D). 
Proof. Let G = (V, E) be a digraph. We have 

Z A ,d(G)= ^ W ^2 II D o(v)a{v) = w - N A,d(G,w). 

w£W a (G) <t£C a (G,w)v&V w&V a (G) 

As the cardinality of Wa(G) is polynomial in the size of G this proves EVAL(j4, D) < 
COUNT(A, D). For the backward direction let G^ denote the digraph obtained from G 
by replacing each edge with t copies of it. We have 

Z AD (G {t) )= ( U A -(nMv)) U D a{v)a(v) = w l ■ N A . D {G,w). 

<j:V->[m] \uv&E / v&V w£W A (G) 

Using an EVAL(j4,D) oracle, we can evaluate this for t = 1, . . . ,\W A (G)\. Therefore the 



values Na,d{G,w) can be recovered in polynomial time by Lemma 4.9 



□ 



4.3.1 The Proof of Lemma [331 

To measure the time needed to perform operations on algebraic numbers, we need to introduce 
their description length. For an integer a G Z let its description length be len(a) := |~log 2 |aj|] + 
2. For a rational number q = | with a, b G Z we define len(g) := len(a) + len(6). For a vector 
a G Q™ we define len(a) = len(ai) + . . . + len(a n ). This defines the description lengths of 
algebraic numbers in Q(0) in standard representation. The description length len(yl) of a 
matrix A G Q nxn is the sum of the description lengths of it entries. 

Let B = {b\, . . . , b n } C Ca be a set of algebraic numbers. By (B) we denote the multi- 
plicative span of B that is (B) = {( • UbeB ftAfc I A 6 G Z, for all b G B, ( is a root of unity} . 
The set B is called multiplicatively independent, if for all Ai, . . . , A n G Z the following holds: 
if b^ 1 ■ ■ ■ b^ n is a root of unity then Ai = . . . = A n = 0. In all other cases we say that B is 
multiplicatively dependent. We say that a set S is effectively representable in terms of B, if for 
given x G S we can compute Ai, . . . , A n G Z such that x ■ b^ 1 ■ ■ ■ b^ n is a root of unity. A set 
B is an effective representation system for a set S, if S is effectively representable in terms of 
B and B is multiplicatively independent. We need a result from [Ri cOlj : 
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Lemma 4.11 (Theorem 2 in [RicOl]). Let a\, . . . ,a n £ Q(9) given in standard represen- 
tation, each of description length at most s. There is a matrix A S ^ nx " such that, for vectors 
A G Z n we have 

n 

Y[ a Ai is a root of unity, if, and only if, A ■ A = 0. (14) 

The description length of A is bounded by a computable function in n and s. 
This theorem has a straightforward algorithmic consequence. 

Corollary 4.12. Let ai,...,a n E Q(9) be given in standard representation. There is an 
algorithm which decides if there is a non-zero vector A = (Ai, . . . , A n ) G Z™ such that 

n 

Y\ a i 4 i s a ro °t of unity. (15) 
i=i 

Furthermore, if it exists, the algorithm computes such a vector A. 

Lemma 4.13. Let B = {bi, . . . ,b n } C Ca be a multiplicatively dependent set of algebraic 
numbers. Then there is a set B such that \B'\ < \B\ and B is effectively representable by B . 



Proof. Since B is multiplicatively dependent, Corollary 4.12 implies that we can compute a 
non-zero vector A G Z n such that b^ 1 ■ ■ ■ b^ n is a root of unity. We can easily make sure that 
at least one of the Aj is larger than zero. Assume therefore w.l.o.g. that Ai > 0. Fix a set 
B' = {b' 2 i • ■ • , b' n } where each b\ is a Ai-th root of bi, that is (&'j) Al = fy. Then 

(n \ A i n-1 

j[l(^) Ai ^ s a ro °* °^ unr ty an d hence b n ■ JJ(^) Al is a root of unity. 
i=2 / i=0 

All operations are computable and effective representation of B by B' follows. n 

Lemma 4.14. Let S C Ca be a set of algebraic numbers. There is an effective representation 
system B for S which can be computed effectively from S. 



Proof. Let B' = S and apply Lemma 4.13 recursively on B. Since the empty set is multi- 
plicatively independent, after at most finitely many steps, we find an effective representation 
system B for S. □ 



Proof (of Lemma 3.1). Let S be the set of non-zero entries of A. By Lemma 4.14 we can 
compute an effective representation system B for S. However, with respect to our model of 
computation we need to be a bit careful, here: assume that S C Q(#) for some primitive 
element 9. The application of Lemma 4.14 does not allow us to stipulate that B C Q(#). 
But in another step of precomputation, we can compute another primitive element 9' for the 
elements of B such that B C Q(9'). Then we may consider all computations as taking place 
in<Q>(0'). 

Assume that B = {bi, . . . , b n }, then every non-zero entry of A has a unique computable 
representation 

n 



Aij = Cij • Y\_ bu iv for Cij a ro °t of unity. 



u=l 
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Let pi, . . . ,pb be f3 = \B\ distinct prime numbers and define A' as the matrix obtained from 
A by replacing in each non-zero entry Aij the powers of b E B by the corresponding powers 
of primes, that is, 

n 



u=l 



Recall the definition of Wa{G) in equation (13). For each w E Wa(G) we can, in polynomial 
time compute a representation w = YlijA^ as powers of elements in S. The effective 
representation of S in terms of B extends to Wa(G) being effectively representable by B. 
Moreover, as S depends only on A, the representation of each w 6 Wa(G) is even polynomial 
time computable. We have 

Z a ,d(G)= Yl w-N a ,d(G,w) 

wdW A {G) 



In particular, for each w E Wa(G), we can compute unique A^i, . . . , \ W;n E Z such that 
w ■ b l 
we have 



w ■ • • • bn w '" is a root of unity. Define functions / and g such that for every w E Wa(G) 



f(w) = Y\_pt w '" and g(w) = b„ w ' v '. 
z/=i i/=i 

Thus we obtain 



^-^ q(w) 

This yields a reduction for EVAL(j4', D) < EVAL(A, D). The other direction follows by 



Z A d(G)= E w ' -jtI-Na'AG^')- 



w'eW A ,(G) 



4.4 Pinning Vertices 



In this section we will prove the Pinning Lemma 3.2 Before we can do this we need to 
establish an important property of partition functions Za,d'- they allow to reconstruct the 
pair (^4, D), in the sense that the values Za,d(G) for all digraphs G determine the pair (^4, D) 
up to isomorphism. 

The notion of isomorphism employed here is given by considering A as (the adjacency 
matrix of) a weighted graph. Let A,D£ £, mxm an d A', D' E C m xm be matrices such that 
D and D' are diagonal. Then A and A' are isomorphic, if m = m! and they admit a bijection 
a : [m] — > [m'] such that for all i, j E [m] we have A^ = A',^ ,.y We call a an isomorphism. 
If further Da = D',^ /a for all i E [m] then the pairs (A, D) and [A' , D') are isomorphic. 

4.4.1 The Reconstruction Lemma 

It will be convenient in the following to fix pinnings and consider digraphs compatible with 
these. To define this, fix <p : [k] — > [m] to denote our pinning. A k-labeled digraph G is then 
a digraph with k distinct vertices labeled 1, . . . , k, in this way (j) is compatible with every 
fc-labeled digraph. The main technical result of this section, is the following. 
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Lemma 4.15 (The Reconstruction Lemma). Let A G C mxm and A' G C m ' xm ' ae imn- 
free, D G jj mxm an( ^ j^i g jgm xm di a go na l matrices of positive vertex weights and <fi : [k] — >■ 
[m] and 0' : [£;] — > [m!\ pinnings. Assume that 

Za,d{4>i G) = ZA' t £j'((j) , G) for all k-labeled digraphs G. 

Then there is an isomorphism a between (A, D) and (A' , D') such that (ft = <j) o a. 

Similar results been given first by Lovasz |Lov67| for the case that A is the adjacency matrix 
of a graph and in [Lov06j for symmetric real-valued A. Schrijver [Sch09j gave a proof for 
symmetric complex-valued matrices. 

We will prove Lemma 4.15 by extending the proof in [Lov06] . We will therefore fix m > m! , 
twin-free matrices A G c mxm and A' G C m ' xra ' and diagonal matrices D G R mxm and 
D' G W 71 xm of positive vertex weights. For convenience we will consider /c-tuples x G [m] k as 
pinnings x : [k] — > [m] and for any x G [m] the expression xx then denotes a (/c + l)-tuple with 
last component x. We define an equivalence relation between tuples x G [m] k and y G [m'] k 
by 

x ~ y iff Za,d(*-, G) = Za',d'{y-, G) f° r au ^-labeled digraphs G. 

The product G\G 2 of two fe-labeled digraphs is the digraph obtained from G\ and G2 by taking 
their disjoint union and identifying vertices with the same label. For A;-labeled digraphs G\ 
and G2 and every x G [m] k we have 

Z A , D (x, GiG 2 ) = Z A , D (x, Gi) • Z A>D (x, G 2 ). (16) 

Lemma 4.16. Let x G [m] k and y G [m'] k such that x ~ y. Then for every and x G [m] 
there is a y G [m'] such that xx ~ yw. 

Proof. Note first that 

m 

Za,.d(x, G) = ^ D^Za^xx, G). 
The condition that x ~ y therefore implies that, for every (k + l)-labeled digraph G 

m m' 

Y, D xx Z A , D (xx, G) = £ D' yy Z A ,, D ,(yy, G). (17) 

x=l y=l 

Define pairs (X\, Yi), . . . , (X s , Y s ) of sets in the following way. The sets X\, . . . , X s partition 
[m] with possibly some empty parts and Y±, . . . ,Y S are partitions of [m'\ with possibly empty 
parts, as well. The partitions satisfy, for every i G [s], that x,x' G X^ if, and only if, 
Za,d(^-x,G) = Za,d(x-x' ,G) for all k + 1-labeled graphs G. The analogue holds for the 
y, y' G Y{ with ZA>,D'(yy,G) = Za\d'{yv' , G) for all k + 1-labeled graphs G. Further for all 
i E [s] , the set Xi U 1^ is non-empty and for all i£lj and y G 1^ we have 

Z^£)(xi, G) = ZA',D'(yy, G) for all fc + 1-labeled digraphs G. (18) 

Define for all i G [s] a (A; + l)-tuple Zj and a function Z(zj,G) such that Zj := for some 
x G and Z(zj,G) := Za,d{ x x,G) if is non-empty and Zj := yy for y £ Yi with 
Z(z.j,G) := ZA',D'{yy-,G) otherwise. Let furthermore 
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Equation (17) now simplifies to 

s 

= Y j c 1 Z(z u G). 

i=i 

We claim that Cj = for all i E [s]. If this holds the proof of the lemma follows as by the 
non-negativity of the vertex weights a = implies Xi ^ and Y{ ^ which, by equation 



(18), provides us with the desired extensions of x and y. We will show, more generally, that 



for all sets I C [s] the equation 

= CiZ(zi, G) for all digraphs G. (19) 



implies that c% = for all i E I. We will prove this by induction on the cardinality of I. 
For I = {i} let Ek+i be the k + 1-labeled digraph on k + 1 vertices which contains no 



edges. If Zj = xi, we have Z(zi,Ek+i) = Za,d{zi, Ek+i) = 1 and thus equation (19) implies 
Ci = 0. An analogous argument proves the case that z, = yy. 

Let |/| > 1. By definition, for all i ^ j E I there is a digraph Gij such that Z(zj,G ! y) 7^ 
Z(zj, Gy). Partition I into equivalence classes Jo, . . . , Jf such that a,b £ J u iff Z(z a , Gy) = 
Z(zb,Gij). Define, for each i E [0, t] the tuple z^ such that z\ = z a for some a E Jj. Let H p 



denote the k + 1-labelled graph product of H with itself taken p times, then by equation ( 16 ) 



and (19) we have 



By definition, the values Z(z' u ,Gij) are pairwise different. It may be the case that one of 
these is zero. Assume w.l.o.g. that = Z(z' ,Gij). For p = l,...,t, this gives rise to a 
homogeneous system of linear equations with a Vandermonde determinant: 

= ^Z(z'„G^)M^c^(z M ,G)). 

f=l \M6J„ / 

Solving this system yields = 2~^eJ,y c Mj Z(z fM ,G) for all 1/ > 1. Since |J^| < n, the induction 
hypothesis implies that we have q = for all i E J v and every u > 1. Therefore, by equation 



(19), we further obtain 

0= £c,-Z(z,-,G). 

which, again by the induction hypothesis, implies c,- = for all j E Jo- n 

Lemma 4.17. Every map tp : [m] — )• [m'] suc/i f/iai Ay- = A^/a^,^ f or a ^ hj £ [ m ] *s 
injective. 

Proof. Let J be the image of -0 and define a mapping t// : J — > [m] such that for each j E J 
we fix an i E [m] with ^(z) = j and define "i/>'(j) = Define <^> : [m] — >• [m] by = '0' o ?/>. 

The precondition of the lemma implies that Aij = Am^m^ for all i,j E [m]. We will 
show first that 4> has a power s which is idempotent, that is (/> s = cj) 2s . To see this, note that 
by the finiteness of [m] not all powers of <j) can be different. Hence, there are i, k such that 
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$ = (j) l+k = (j) l+ik for all i G N. Multiplying from the left with <jj we have 4> l+j = 4> l+ i +lk 
for all i,j G N. Let i be large enough such that ik — I > then with j = ik — I we have 
2(Z + i) = I + j + ik. Hence = <^+-»'+ ifc = (f> 2( - l +ti . 

Now j4jj = Afpu^fj) implies Aij = A^u^su^ for all s. Hence for s such that 4> s = 4> 2s we 
have 

= ^(t> s {i)<t> s {j) = ^V 2s (i)<£ s 0') = ^4> s (i)j ^ or a11 G [ m l- 

Recall that ^4 is twin-free. Thus (j) s is the identity and therefore (f> must be bijective. Further 
by the definition of this implies that if) is injective. n 

Lemma 4.18. Let x G [m] m and y G [m'] m such that x is bijective and x ~ y. TTjera i/iere 
is an isomorphism a between (A, D) and (A', D') such that y = a o x. 

Proof. Define Cry as the m- labeled digraph on m vertices with the single edge ij. Assume 
w.l.o.g. that Xi = i for all i G [m]. By x ~ y we have 

= Z A>D (x, dj) = Z A > tD >(y, Gij) = A'y.y.. (20) 



Since A' is twin- free Lemma 4.17 implies that y is injective and by the condition that m > m! 
it is thus bijective and m = ml . This gives rise to an isomorphism a : [m] —> [m] between A 
and A' by defining a(xi) = a(i) =: yi for all i G [m]. 

It remains to show that Da = D',^ holds for all i G [m]. Let Q = {i G [m] \ A^* ^ 0} 
and P = {j G [m] \ A* t j ^ 0} be the set of indices corresponding to non-zero rows and 
columns of A. If one of these sets is empty, then Q = P = and A = (0) since A is twin-free 
and then the following claim implies D\\ = D'^. 

Claim 1. We have 

m m 

i=i j=i 

Proof. Let Go be the m-labeled digraph on vertex set [m]U{u} without edges. By x ~ y we 
have 

m m 

E" ^ = Z A;D (x,G ) = Z A > tD i(y,G ) = E D ' a (j)a(j)- 

i=i i=i h 



Assume therefore that Q,P 7^ 0. Let p : Q — > N be an mapping, we say that it is non- 
vanishing if there is an i G Q such that p(i) > 0. Define G^ as the m-labeled digraph on 
vertex set [m]0{v} with edges E = {{iv) p ^ \ i G Q}, that is for each i G Q there is an edge 
with multiplicity p(i). Note in particular that the edge iv does not exist if p(i) = 0. 
For every non-vanishing p : Q — > N, we have 

Z A>D (x,GW) = j2D jj H(A ij r. 

jeP i&Q 

By x ~ y we have Z a ,d(x-i G^) = Z A / t jy(y, G^ 1 ') and therefore, 

e D n ri(^) p(i) = e n«(i)«(i)) pW 

jeP ieQ jeP i&Q 
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Now a is an isomorphism between A and A', that is A^ = ^(iWi) anc ^ thus the above 
simplifies to 

= E [Pii ~ D kjMJ)) II ( Ai ^ Pil) for a11 non- vanishing p:Q^N. (21) 
We shall show the following 

Claim 2. Let JCPof cardinality at least 2 and assume that 

= E { D n ~ D '«UMj)) H( A ij) P(t) for a11 non-vanishing p : Q -> N. (22) 

Then for every j G J, there is a set J' C J with j G J' such that 

= E ( D n ~ D *U)a(j)) II (^i) P(i) for a11 non-vanishing p : Q -> N. (23) 
jeJ' ieQ 

Before we prove this claim, let us see how it finishes the proof of the lemma. Its statement 
will be particularly true for J = P and the condition of which is given by equation (21). 
Thus, for every j £ P iterative application of this claim eventually implies 

= (Djj - D' a(j)aU) ^ YliAijf^ for all n on- vanishing p : Q N. 

ieQ 

Since j G P there is an i G Q such that Aij ^ 0. Fix a p : Q — > N which satisfies p(i) = 1 and 
p(i') = for i' ^ i Then we have = (^Djj - D^. )a ^j Ay and thus £>-,-.,• = D' a ^ )a{j) . 

Therefore we obtain D,-, = P/ , s , .n for all j G P. If A contains no zero columns, we 
are done as then P = [m\. Otherwise, since A is twin-free, there is exactly one k such that 
^*,k = and PU {A;} = [m\. Thus Claim jlj further implies P/^ = D' a ^ a ^ k y 

The Proof of Claim [2| As J is of cardinality at least 2, we may fix distinct G J. By 
A being twin-free we may fix an i* G Q such that Ai*j ^ A^y. Let p : Q — > N be a non- 
vanishing mapping. Define a non- vanishing p q : Q — > N by p q (i*) = p(i*) + q and = p(i) 
for all i ^ i* . The condition of the claim implies that 

o = E(^-^o>o))ri(^) P9W 

jeJ ieQ 

= E(^) 9 (^i-^(iw,))II(^i) Ki) 

jGJ ieQ 

Partition J into equivalence classes Jo, . . . , Jt such that Ai*j = Ai*ji iff j,j' G Jk for some 
A; G [i]. Define, for all k G [0,t] values at '■= A{* u for some v G J^. Then 



k=o \jeJ k ieQ 
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It may be the case that one of the a& values is zero. Assume w.l.o.g. that ao = 0. By condition, 
the above holds for all q. In particular, for q = 1, . . . , t this gives rise to a homogeneous system 
of linear equations 



i 



= Ew1E^r% (J ))n^) p( 

k=l \jeJ k i<=Q 

Since all the values a±, . . . are pairwise distinct and non-zero, this system has a Vander- 
monde determinant and is thus solvable, which implies, for all k 6 [t], 

0= Y1 { D 3J - D a(j)a(j)) I[( A *j) Pii) for a11 non-vanishing p : Q -> N. 
In combination with the condition ( |22[ ) of the claim we also obtain 

= Yl { D rJ ~ ^aOXi)) II(^) P(i) for a11 non-vanishing p : Q -> N. 

Since the Jo , . . . , Jt are a partition of J and each of them has cardinality strictly smaller than 
J, the claim follows. n 

Now we are able to prove the Reconstruction Lemma. 



Proof (of Lemma 4-15). Let A, D be m x to matrices and A',D' be m! x ml matrices as 
given in the lemma and (ft, (j)' pinnings. In accordance with the notation of this section we 
define x := and y := <p' and we assume w.l.o.g. that m > m' . The precondition that 
Za,d(4>i G) = Za',d'(4>' ') G) for all digraphs G is tantamount to x ~ y. 

It will be easier to work with surjective pinnings. Extend x to a surjective x' := xz € [m] 1 



with / > m. Lemma 4.16 implies the existence of a z' such that for y' = yz' we have x' ~ y'. 
Observe that if we find an isomorphism a between (^4, D) and (A' , D') such that y' = aox' 
then the proof follows as this implies y = a o x. 

We shall prove the existence of such an isomorphism. Assume w.l.o.g. that x\ = i for all 
i £ [m] — we can achieve this by permuting the labels appropriately. Define x'|/ = (a^)jgj 
for a set of indices / C [I]. 

Claim 1. If x"x ~ y"y for some x £ [m] and y £ [m'\ then x" ~ y". 

Proof. To see that this holds, assume the contrary. That is, for some fc-labeled digraph G we 
have Za,d{^' i G) ^ Za' : d'(y" , G). Let G' be the digraph obtained from G by adding a single 
isolated vertex labeled k + 1 , then 

Z a ,d(x"x,G') = Z a ,d(x",G) + Z A >, D '(y",G) = Z A ',d> (y"y, G') 

contradicting the assumption. -I 



4.18 



This claim implies that x'|[ m | ~ y'| [ m ] • The m-tuple x'| r m i is bijective by definition. Lemma 
thus implies y'|[ m ] = a o x'|[ m ] for some isomorphism a between (^4, D) and (A', D'). 

It remains to show that for all j S [m + 1, 1] we have y'j = a(x'j). The bijectivity of x'|[ m ] 
implies x'j = x' r for an appropriate r £ [m\. Define I = {1, . . . ,r — l,r + l, . . . , to, j}, then x'|/ 



is bijective. Therefore, Lemma 4.18 implies the existence of an isomorphism ft between (^4, D) 
and (A',D') such that y'|/ = /3ox'|/, particularly implying bijectivity of y'|/. It follows thus 
for our original automorphism a that y'j = y' r = a(x' r ) = a{x'-). n 
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4.4.2 The Proof of The Pinning Lemma. 



Applying the Reconstruction Lemma 4.15, we will now prove the Pinning Lemma. 

Lemma 4.19. Let A G C™ xm be a Hermitian matrix and D G M™ xm a diagonal matrix of 
positive vertex weights. Let [A] be the twin resolvent of A. Then 

EVAL pin ([A],D^) = EVALP''"(A, D) and EVAL(L4], D^) = EVAL( J 4, D). 



Proof. Let r be the twin resolution mapping on A. By the Twin Reduction Lemma 4.5 2) 
we have 

Za,d{4>-, G) = Zrju D [A\ (r o 4>,G), for every digraph G and pinning <f>. 

This yields EVAL pin (A, D) < EVAL pin (L4], D^). On the other hand, for some instance G 
and (f)' of EVAL pin (L4], D^) let 4> be a pinning on A such that to cfy = <f> , then the above also 
yields the backward direction. The result for EVAL(j4, D) follows as the above holds likewise 
for empty pinnings. n 

Lemma 4.20. Let A G C™ xm be a Hermitian and twin-free matrix and D G M™ xm o diago- 
nal matrix of positive vertex weights. Then 

EVAL p/ "(A D) = EVAL(A D). 

Proof. It suffices to show EVAL pin (^, D) < EVAL(A, D) as the other direction holds trivially. 

Let G = (V, E) and a pinning </> be an instance of EVAL pm (^4, D). By appropriate permu- 
tation of the rows/columns of A and D (cf. Lemma 4.3) we may assume that [k] = img C [m] 
for some k < m. Let G = (V, E) be the graph obtained from G by collapsing the the sets 
for all i £ [k]. Formally, define a map 



7(f) 



i ,v G (j) for some i G [k] 
v , otherwise 



Then G is a multi-digraph (with possibly some self-loops) defined by 

V = [k] U (V \ def(0)) 
E 1 = {7(14)7(2;) | lif G -B}. 

Hence (</>, G) = ^4,£> (idpy ; G) where idj fc ] denotes the identity map on [k]. Call two 
mappings x> "0 : [&] — * [ m ] equivalent if there is an automorphism a of (^4, D) such that 
X = a. o ■0. Partition the mappings t/> : [A;] — > [m] into equivalence classes Ji, . . . , I c according 
to this definition and for all i G [c] fix some ipi G Li. Assume furthermore, that ip± = idny. 
Clearly for any two Xi^ from the same equivalence class, we have Za,d{XiF) = Za,d{iP, E) 
for every graph F. Therefore, for every graph G' , 

= £ Z A , D (^, (?') • f J2 II ^WW I ( 24 ) 
i=i VV'S/i ^edcf(V') / 



Define, for each z G [c] the value c; = ( X^eA ELedef(i/>) A/>0)V>(t>) J • We claim the following 
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Claim 1. Let I C [c] be a set of cardinality at least 2 such that 1 G /. Assume that we can, 
for every ^-labeled digraph G", compute the value 

Y,ci-Z A>D ^ U G'). (25) 

Then there is a proper subset I' C I which contains 1 such that we can compute, for every 
A;-labeled digraph G" , the value 

Y,Ci-Z A)D {^,G"). (26) 
iei' 



This claim will allow us to finish the proof. To see this, note first that by equation ([24]) we 
can compute the value ( |25[) for / = [c] and G' = G. Thus after at most c iterations of Claim [l] 
we arrive at c\ ■ Za,d{iPi,G). Further, c\ is effectively computable in time depending only on 
D and therefore we can compute Z a> d (idru , G) = Za,b{4>^G). This proves the reducibility 
EVAL pin (^4, D) < ENAL(A,D). 



Proof Of Claim [TJ Assume that we can compute the value given in (25). Lemma 4.15 



implies that for every pair i ^= j E I there is a ^-labeled digraph T such that 

Z A)D (tl> it T) ^Z A7 D&j,T). (27) 

Fix such a pair i 7^ j G I and a graph T satisfying this equation. Note that this graph can be 
computed effectively in time depending only on A, D and tpi,tpj- Let G s denote the digraph 
obtained from G by iterating s times the fc-labeled product of G with itself. We can thus 
compute 

c 

a ■ z a ,d(^ g't s ) = ciZ AtD &i, g') ■ z AiD (^,ry. (28) 

Partition / into classes Jq,...,J z such that for every v G [0, z] we have G J v if, and 
only if, Za,d(}1>v >T) = Za,d (Vy ; F) • Since one of these sets J y contains 1 and all of these are 
proper subsets of /, it remains to show that we can compute, for each v G [0, z], the value 



ciZ A ,D{^i',G'). 



To prove this, define x u := Za _d(V , « / ) T) for each v G [2] and an i' G J„. Equation (28) implies 
that we can compute 

f / 

Ci'Z A ,D(ipi>,G') 

One of the values x v might be zero. Assume therefore w.l.o.g. that xo = 0, then evaluating 
the above for s = 1, . . . , z yields a system of linear equations, which by Lemma |4.9| can be 
solved in polynomial time such that we can recover the values ^i'eJ„ c *' 'Z A ,d{^v j G') for each 



v>\. Using equation (25) we can thus also compute the value 



^Cj- Z AiD (ipi,G') ) X! c i' z A,D(^i',G') = ^ Ci>Z A}D (i/ji,,G'). 



4.5 Basic Complexity Results for Congruential Partition Functions 
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Proof (of t he Pi nning Lemma 3.2). Fix the twin resolvent A' 
By Lemma 



4.20 



we have EVAL pm (A', D') = EVAL(A', D') 



[A] of A and let D' = D^ A l 
Hence together with Lemma 



4.19 



we obtain the chain of reductions 

EVAL pin (A D) = EVAL pin (A', D') = EVAL(A', D') = EVAL(A, D). 



4.5 Basic Complexity Results for Congruential Partition Functions 

A Basic Polynomial Time Case. We provide the following basic tractability result which 
is a straightforward extension of Theorem 6 in [BG05J. 

Lemma 4.21. Let A £ (rynxm, ^ g a }j erm iti an matrix and D a family of diagonal matrices 
of vertex- weights. If A has rank 1 then EVAL P '"(A, D) is polynomial time computable. 



Proof. Let G = (V, E) be a given digraph and 4> a pinning. Let V = V \ def(</>) then 

[8(v)J 
r(v),cr(v) 



z A ^ G )= £ n^),wii<;!( 



rpC(T:V^[m] uv£E vGV 



As rank A = 1, there are vectors a, b E such that A = ab T . Then, for every configuration 
a : V — > [m], 

nA TT TT de § + ( v ) u de s~ («) 

Aa( u ), a ( v ) =11 a a(u) b a{v) = [[ a a(v) b a(v) . 



Therefore, 



Z A>D ^G) = Y, Il^HllCw 

4>Ca:V^[m] uv&E v&V 



TT ^ {v) b d : t \ {v) T TT a de ^ (v) b d % {v) D ld ^\ , 

11 4>(V) <p( v ) I / J 11 V\ v ) <J \ V ) <T(V),<T(V) 

vedei((j)) ) <r.V'^[m]veV 

m 

de g + ( v ) 5 de S~ ( v ) 
rj>(v) <f>(v) 
«Sdef(0) veV i=l 



n deg+(^),deg (v) TT deg+(i;),deg («) n {d(v 

a <j>(v) °<t>{v) 11 Z^ a i °i 
Edef(0) veV i=l 

And the term in the last line can be evaluated in polynomial time. 



Basic #P-hardness Results. We need the following extension of the ^tP-hardness crite- 
rion of Theorem 11.11 

Lemma 4.22. Let A € M™ xm be a non-negative symmetric matrix and D € IR™ xm a diagonal 
matrix of positive vertex weights. If A contains a block of rank at least 2 then EVAL(A,D) 
is #P-hard. 

We will give the proof of this Lemma below. From this result we derive a basic #P-hardness 
criterion for congruential partition functions. 

Lemma 4.23. Let A be a Hermitian matrix and D = (-D^) c ez w such that I)M is a diagonal 
matrix of positive vertex weights. If abs (^4) contains a block of row rank at least 2 then 
E\AL pir '(j4,D) is #P-hard. 
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Proof. Note first, that we have EVAL pin (abs (A) {2) , £>I°J) < EVAL pin (A, D). To see this, 
let a graph G = (V, E) and a pinning 4> be an instance of EVAL pin (abs (A), flM). Let 
G' be the graph obtained from G by adding for each existing arc uv an arc uu, we have 

As abs (A) c ontains a block of row rank at least 2 the matrix abs (A) ^ does so as well. 



By Lemma 4.22 the problem EVAL pin (abs (^) (2) , £> M ) is #P-hard. 



The proof of Lemma 4.22 The proof follows from Theorem 1 1 . 1 1 and the following Lemma 



which is a straightforward extension of Theorem 3.2 in [DGOO]: 



Lemma 4.24. Let A G M.™ m be a symmetric matrix with non-negative entries such that 
every pair of rows in A is linearly independent. Let D G M™ xm be a diagonal matrix of 
positive vertex weights. Then 

EVAL(A) < EVAL( J 4, D). 

The proof adapts the one in |DGO 0J : we present it for completeness. The following lemmas 
are restatements of those in [DGOOJ (see Lemma 3.4, 3.6, 3.7 and Theorem 3.1). 

Lemma 4.25. Let A G M™ xm be symmetric and non-singular, G = (V, E) a graph and 
FOE. If we know the values 

f r (G) = £ c A (a) J] Al {uMv) (29) 

cr:V^[m] uv£F 

2 

for all r G [(\F\ + l) m ], where ca is a function depending on A but not on r. Then we can 
evaluate 

° A ( a ) II ^m)a{u)a{v) (30) 
a:V— >[m] uv£F 

in polynomial time. 

Proof. As A is symmetric and non-singular, there is an orthogonal matrix P such that 
P T AP =: D is a diagonal matrix with non-zero diagonal. Every entry of A r satisfies 
Alj = {PD r P T ) i3 = £™ =1 P a{u)fl P a{v)tl (D^y. By equation §j} 



f r (G) = Y, C ^ IT T, P °(u)H P <r(v)»(D^y 
a:V^[m] UV&F /Lt=l 

Define W = {YYiLi(Dii) 011 I < on for all i G [m], Y^=x a i = 1-^1} which can be constructed 
in polynomial time. We rewrite 

fr(G) = CwW r - 

wew 



for unknown coefficients c w . By interpolation (cf. Lemma 4.9), we can recover these coeffi 



cients in polynomial time and can thus calculate fo(G) = Su»eW c ^- We nave 

m 

f (G) = Yl C ^°) II J2 P ^Pa(v)^D^)° 

<r:V->[m] uvEF H=l 



which is equal to ( 30 ) by inspection 



4.5 Basic Complexity Results for Congruential Partition Functions 
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Lemma 4.26. Let A £ ^ mxm £ e a symmetric matrix in which every pair of distinct rows 
is linearly independent. Let D £ jj mxm ^ e fl d{ a g 0na i matrix of non-negative vertex weights. 
Then every pair of rows in ADA is linearly independent. Furthermore there is an < e < 1 
such that 

\(ADA)ij\ < e^J (ADA)u(ADA)jj 

Proof. Define Q = AD^I 2 \ We have ADA = AD^^D^ l /^A T = QQ T . That is (ADA)ij = 
{Qi,*i Qj,*) every pair of rows in Q is linearly independent as it is linearly independent in A. By 
Cauchy-Schwarz (Qj,*, Qj,*} < \j (Qi,*, Qi,*)(Qj,*, Qj,*}, which implies that the corresponding 
2x2 submatrix of ADA defined by i and j has non-zero determinant. The existence of e 
follows. n 

Lemma 4.27. Let A £ l mxm be a symmetric non-negative matrix in which every pair of 
distinct rows is linearly independent. Let D £ ]j mxm fr e a diagonal matrix of positive vertex 
weights. Then there is opGN such that the matrix (ADA)^ is non-singular. 

Proof. Let A' = ADA and consider the determinant 

m 

det(A')= £ ±JI 

TTSSm »=1 

Let S m be the set of permutations it : [m] — > [to]. For some ir £ S m define t(ir) = \{i \ ir(i) ^ 
i}\. Let e be as in Lemma 4.26 Then 



■m 



i=l i=l i=l i=l 

Let id denote the trivial permutation. Then 

(m \ P / m s 

«=1 / 7reS m \{id} \i=l / 



By equation (31), we have 




and hence, as < e < 1, the matrix (ADA)^ is non-singular for large enough p. 



4.27 



Proof (of Lemma 4-^4)- Fix a a p £ N such that (ADA)^ is non-singular by Lemma 
Let a graph G = (V, E) be an instance of EVAL(A). Define a graph G' = (V',E') by 
V = {u , . . . ,u d -i I u G V,d = d G (v)} and £' = E" U for sets E", E'" to be defined 
next. Let E'" be the set which contains, for each v £ V the cycle on {vq, . . . ,Vd G ( v )-i}, i.e. 
{vo^ij • • • j v d G (v)-i v o} E'" ■ Let E" be the set of edges, such that each edge from E incident 
with v is connected to exactly one of the «j. 

We further construct a graph G' p r from G' as follows. Let T p be a graph which consists 
p many length 2 paths connecting terminals a and b. Let T Pir be the series composition of 
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r many copies of T p . G' pr is obtained from G' by replacing each edge in E'" with a distinct 
copy of Tp <r . We call a and b the "start" and "end" vertex of T p and T p<r has start and end 
vertices induced by the start and end vertices of of the series composition of T p . 

For a graph H with designated "start" and "end" vertex let Z Ai D(i,j',H) := Za,d(,4>iH) 
such that (ft pins the start vertex to i and the end vertex to j. 

Claim 1. Let C = (ADA)^\ then with X = D^/ 2 \ we have for all i,j 6 [m] and r 6 N, 

Z AD (i,j;T p>r ) = (XuX^-^XCXYij. (32) 

Proof. Straightforwardly, 

Z A>D (i,j;T p ) = \Y^A lk A kj D kk \ = (ADA)$ = dj 

and therefore 

r r 

ZA,D(i,j;T p , r ) = ^2 Y\Z AiD (a(k),a(k + l);T p )Y[D a{k)ia(k) 

cr:[r+l]-+[m] k=l k=2 
<r(l)=i, <r(r+l)=j 

r 

= (XuXjj) 1 ^2 ]^[ Xa(k),a(k)C a (k),a(k+l)Xa(k+l),a(k+l) 

<r.[r+l]-y[m] fc=l 
tx(l)=i, tr(r+l)=j 



By inspection, the last line equals the right hand side of equation (32) — as claimed. 



For convenience, we count the indices of the vertices i?o, . . . ,Vd G ( v )-i modulo dc(v). Particu- 
larly, Vd G ( v \ = vo and we have, for all r G N, 

d G (v)-l 

z A,D(G'p, r ) = XT II A a(u)(r(v) IJ At(u),<t(») JJ JI Z A)D (a(vi), o-(v i+1 ); T p>r ). 

ar:V'-+[m]uveE" veV ueV i=0 

As the vertices in 1/' are grouped according to the vertices in V we have 

d a (v)-l 

n = n n a^m^- 



By Q, the expression ILeV' Ar(«VM TLev E[f=o 1 -^4,d(o"(^), c(^+i); r P)J .) equals 



d G («)-i (xcxy d G (v)-i 

n n b %w««)v ; : - v 3 ] — ; = n n (^^^w^i) 

v£ V i=0 yV o-(wi) ) ff(«i) J<1 -o-(«<+i),<T(o»+i) „ e y i=o 



Thus 



^A,D(Gp >r ) = n ^(«)^«) n ( xc ' x )a(«)a( 1 ,) 

cr:V'-+\m\uveE" uveE'" 
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Given the EVAL(t4, D) oracle, we can, in polynomial time, evaluate this expression for every 



r which is polynomial in the size of G. Lemma 4.25 implies that we can compute the value 



z - yi n 4t(«mio n ^(uMv) 

o-.V'->[m] uv&E" uv&E'" 

It remains to show that Z = Za{G). To see this, note that, for every configuration a : V — > 
[m] the weight n^es" At(m)o-(u) Yiuv&E 1 " ^a(u)a(v) i n the above expression is zero unless the 
following holds: For all v E V we have a(vo) = . . . = a{vd-i) for d = do{v). For such a 
configuration, define a' : V —> [m] by o~(v) = ct(vq) for every v £ V. Then 

uv<EE" uvEE 

Since every configuration a' : V — > [m] arises this way this finishes the proof. n 

5 Connected Hermitian Matrices 



In this section we will prove Lemma 3.4 To do this we will split the proof into two parts 
- one for the non-bipartite case and one for the bipartite case. The following two lemmas 
provide these two parts. 

Lemma 5.1 (The Non-Bipartite Case). Let A be a connected non-bipartite u-algebraic 
and Hermitian matrix and D a diagonal matrix of positive vertex weights. Then either the 
problem EVAL p '"(j4, D) is jfP-hard or the following holds. 

There is an u-algebraic Hermitian matrix H and a family U = (t/" c ")c6Z„ of diagonal 
matrices which define an (H— STD) problem such that 

EVAL pin (A, D) = EYAL pin (H, U). 

Lemma 5.2 (The Bipartite Case). Let A be a connected bipartite u-algebraic Hermitian 
matrix and D a diagonal matrix of positive vertex weights. Then either EVAL P '"(^4, D) is 
jfP-hard or the following holds. 

There is an ui-algebraic Hermitian matrix A' and a family U = (J7^) c eZ w of diagonal 
matrices which define a (B H STD) problem such that 

EVAL pin {A,D) = EVAL pin (A' , U). 



The proof of Lemma 3.4 follows directly from these two lemmas. In the remainder of this 



section we will give the proof of the Non-Bipartite Case Lemma 5.1 The proof of the Bipartite 



Case Lemma 5.2 is conceptually very similar, therefore we omit it, the full proof can be found 
in |Thu09j . 



The proof of Lemma 5.1, The strategy is to filter the problems EVAL pm (A, D) until we 
are left with only those for which we can find a polynomial time equivalent (H— STD) problem 
EVAL pm (ff, U). This process of filtering will be performed by transforming in several steps 
the problem EVAL pin (A, D) into a problem EVAL pin (C, D) where the matrix C and the family 
D satisfy certain conditions, which will be given below. In each step, if a transformation is 
not possible, it will be shown that this implies #P-hardness of the corresponding problem. 
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A bit of preparation is necessary. We say that a family D = (£)M) c6 z w fits a matrix C if 
the dimensions of each £H C J match those of C. The matrix C under consideration will have a 
tensor product structure X (g) H. Recall that we index the entries of such a matrix by pairs 
(fj,,i),(v,j) such that C(^,i),(wj) = -X^i/ • Hij. The diagonal matrices (Z^W ) ce ^ fit C and 
therefore we index these in the same way. For convenience, we will alter the notation of the 
indices of these diagonal matrices a bit and let 

Shape Conditions. Let C be an w-algebraic Hermitian matrix and D = (£M c J) cg z u a 
family of diagonal matrices of vertex weights that fits C . We define conditions on the shape 
of C and D, two of which determine the shape of C. 

(CI) There are r, m G N, q G {—1,1}, a non-singular normalized matrix 77 G U^ xr and 
positive vectors v,w £ M.™ such that t>i < . . . < t; m and w depends linearly on u and 



q ■ viWiH ... ? • viw m H 



C = q ■ vw T <g> H 



q ■ v m wiH ... q ■ v m w m H 
The submatrices CV^,*)(i/,*) = ? • v^WyH are the £iZes of C. 
(C2) The matrix H is a normalized Hermitian complex Hadamard matrix. 
Three conditions give structure of the diagonal matrices in D. 



(Dl) Z)M has positive diagonal and for all c G Z w we have £>I- C J = £)W. 

(D2) There is an m x m diagonal matrix A^ with positive diagonal such that £)M = 

(D3) For each c G there is a diagonal matrix £/I c l G ({0} UU) rxr and a diagonal matrix 
AW such that DW = A' e ' ® £/M. Further, if flW = then £/H = 0. 

These conditions will play a central role throughout the whole section. Once we make sure 
that they are satisfied, the problem EVAL pm (C, D) has all properties necessary to find a 
polynomial time equivalent (H STD) problem. Then the following lemma provides us with 
the last part of the proof of Lemma |5.1| 



Lemma 5.3. LetC,TJ satisfy conditions (C1) ; (C2) and (D1)-(D3). Let a matrix H and a 
family U = (C/I c H) ce z„ be defined as in conditions (C1),(C2) and (D1)-(D3), respectively. 
Then 

EVAL^CD) = EVAL p/n (#,U). 



Before we give the proof of the Lemma, let us first develop a small technical result which will 
play an important role in that proof. 
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Lemma 5.4. Let A G C mxm and A' G <C m ' xm ' be u- algebraic. Let~D = (£> Ic] ) ceZu &e o/amtli/ 
ofmxm diagonal matrices and D' = (D'I c ]l) ce x„ a family of m! x m' diagonal matrices. Let 
A" = A A' and D" = (ZH C 1 ® D'W ) c ez„ • There are projections r and r' smc/i i/iai /or every 
digraph G and pinning (ft we have 

Za",D"(((>,G) = Z A ,n(TO(f), G) ■ Z A /^'(t' o (p,G). 

Proof. For convenience, we write the indices of A" as pairs according to the tensor product 
structure. That is 

A'(i,k),(j,l) = Aij ' f° r an h3 G { m \ an d k,l G [m'\. 

This holds analogously for the elements of D". We define r : [mm'] — > [m] and r' : [mm 1 ] — > 
[m 1 ] by r(i, k) = i and r'(i, k) = k for all (i, A:) G [m] x [m']. 

Let G = (V, E) be a digraph and <ft a pinning w.r.t. A". Define V = V \ def((/>). 

m*. g) = En n ® ^ ^w^) 

</>Co-:V-»[m]x[m'] uuefi vSV 

y ] ^rocr(u),roo-(u)^4 T 'o(T(u),T'oo-(ii) 
i/iC(r:V->[m]x[m'] uveE 

n D l9{v)i . 
TOff(i)),TOff(t)) r' oo(v),t' oa(y) 

For each configuration a ^ (ft we can form two independent maps a* and cr' such that r o C 
a* = t o a and r' o C a' = t' o o~. Therefore we obtain 

\to0C<t*:V--)-[to] uv&E 1)67' 

e n n ^J^M 

^'o^Cff':F->[m'] u»6E DgV 

In other words, Za",d"((P, G) = Za,l>(t o cp,G) • Za>,y>>{t' o (f>, G). n 
Proof (of Lemma 5.3). Define M = q ■ vw T and let A = (Al c H) cg ^ be a family of m x m 



diagonal matrices defined as in conditions (D2) and (D3). In particular, C = M®H a nd fo r 



each c G Z w we have D^l = A' c ' ® [7^. Then the following claim follows from Lemma 5.4 



Claim 1. There are mappings r and r' such that for every digraph G = (V,E) and every 
pinning 0, we have 

Zc,d(& GO = Z ma (t o<f>,G)- Z h>xj (t' o&G). 



By Lemma 4.21 the function EVAL pm (M, A) is polynomial time computable and therefore 
Claim [l] implies the existence of a reduction witnessing EVAL pin (C, D) < EVAL pin (#, U). 

It remains to prove the reducibility EVAL pin (H, U) < EVAL pin (C, D), the proof of which is 
slightly more complicated. The complications arise from the fact that there may be graphs 
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G for which Zm,a{t </>, G) = implying Zc,o(4>, G) = although Zh,u{t' o <f>, G) might be 
non-zero. 

We will tackle this problem by (p ■ ui + 1) -thickening. Let G^ puj+1 ^ be the graph obtained 
from G by this operation. Note that (p ■ uj + l)d(u) = 9(f) (mod w) for all vertices v whence 
it thus follows by inspection that 

Z h , v (t' o 0, = Z h(p .„ +1)jV (t' o<j>,G). 

Since every entry of H is an w-th root of unity, we have = H. Therefore, 

ZcM G(P " +1) ) = z m,a(t ° 0, G^ +1 >) ■ Z HiV (t' o 0, G). 
We will show the following claim. 

Claim 2. There is a computable pGN which depends only on M and A such that Zh,u (t' o 
(/), G) ^ implies Z A /,a( t o 0, ^ 0. 

Before we give the proof let us see how this finishes the proof of the lemma. Given this p 
we can compute Zm,a.{t <fi, in polynomial time. If it is zero then we know that 

Zh,u( t ' 4>i G) is and we are done. Otherwise, we can compute 

Zma(j ° & G^+^Y 1 ■ ZcM& G (pam) ) = Zh,u(t o</>,G). 



This yields EVXL pin (H, U) < EVAL pin (C, D). 

Proof of Claim|2[ Recall that Z M ,j\( TO( t>, G ) = -^ TO T >A (ro0, G). We will therefore 
focus on Z vw t a( TO( / ) ) G). Define V = V\def(ro0) and recall that by (CI) we have w = A- v 
for some positive A. Therefore 

Z vw T A (ro^,G) = II^)II A K.) 

ro</>Ccr:V->[m] ab£E a£V 

= aw y rr TT v d [ a \A^f, , 

L 4 11 c(o) 11 cr(o) o-(a)<r(a) 

TO<£C<x:V->-[ra] aedef(</>) aeV 

- Al £ l TT v d{a) V TT /^a^J 

~ 11 V TO<t>{a) 11 y <r(a) <x(a)<r(a) 



aedof(0) cr:V'^>H aeV" 



It thus follows that 



w^ 15 ) = Als| n ^jr G(o) e n *) +1HG(a)A £S SG(a)1 

agdef(</>) cr:V'-+[m]a£V' 

m 

= A |jB| Y\ w (pw+1) ' dG(a) ^(ph)+i)-d G (o)^ia G (a)] 

aedef(</>) aSV i=l 

By this calculation, a (t o 0, G^ pbJ+1 ^) is non-zero for some p € N, if for each a € V, 

m 



5.1 Satisfying Shape Conditions (CI) and (Dl) 
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To finish the proof it suffices to prove, for each a G V , the existence of a p a such that with 



P > Pa the inequality in (33) is satisfied. If further each p a depends only on M and A, then 



with p being the maximum of all these p a the proof of the claim follows. 

The assumption Zh,u(t' o 0, G) ^ implies that for each a G V the matrix JJ^ a ^ a ^ ^ 
and thus by condition (D3) the matrix dPgwI contains a non-zero diagonal entry. By the 
General Principles of Section [4. 1| we may assume that G is connected. Further the reduction 
is trivial if G contains no edge, therefore we may assume further that dc(a) > 1 for all a G V. 

Fix an a G V' and let i a be the maximum index i G [m] such that A|f G< - a ^ ^ 0. If there 
is only one index i G [m] such that a|? g ^-" ^ 0, then the inequality in (33) holds for all 
P > Pa '■= 0. Otherwise, let z G [m], with z 7^ i a be maximal such that A|f G ^ a ^ 7^ 0. The 



inequality (33) is satisfied if 



By the definition of z this is the case if 

. \ (pu)+l)-d G (a) 



> 



(p£j+l)-d G (a) a [9o(a 



(34) 



> 



A^ G(a) 



(35) 



As dc{a) > 1, w > 1 and < u z < V{ a , there is a p a such that for all p > p a the inequality 



(35) (and hence (33)) is satisfied. To finish the proof, observe that p a depends only on A and 
vi,...,v m . n 



5.1 Satisfying Shape Conditions (CI) and (Dl) 

We will now show, in a sequence of several small steps, how to transform the initial problem 
EVAL pm (v4, D) into a problem satisfying the shape conditions. In each step, if transformation 
is not possible we will argue that this is due to the problem at hand being #P-hard. 

The first step will provide us with the transformation of EVAL pm (A, D) into a problem 
satisfying (CI) and (Dl). 

Lemma 5.5. Let A G C^*™ he a Hermitian uj-algebraic matrix which is connected and 
non-bipartite and D G W r ^ xn a diagonal matrix of positive vertex weights. Then either 
EVAL pin (A,D) is #P-hard or the following holds. 

There is a matrix C satisfying condition (CI) and a family D of diagonal matrices which 
satisfies condition (Dl) such that 

EYAL pin (A,D) = EVAL p/n (C,D). 

Further, there are two technical results following from the proof. 

(1) H(D^)H T = tr(I>M^)J r for all \i G [m]. 

(2) There are constants di,...,d m such that Z}I°1^ = d^D^' 1 for all 



We split the proof into two parts. The first of which is given in the following lemma. 
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Lemma 5.6. Let A E C^ xn be a Hermitian u-algebraic matrix which is connected and non- 
bipartite and D E R^*" a diagonal matrix of positive vertex weights. Either EVAL pm (A,D) 
is jfP-hard or the following holds. 

There is a Hermitian uo-algebraic and twin-free matrix C and a family D = (-D" c ")ceZ„ °f 
diagonal matrices which satisfies condition (Dl) such that 

EVAL pin (A,D) = EVAL p/n (C,D). 

The matrix C further has the following properties. There is a q E { — 1, 1} and positive vectors 
v, w E such that v\ < . . . < v m and w = Xv for some positive A such that 

I V!WlP 11 ... V X W m P Xm 

C= : •.. : 

\ v mWl P ml . . . v m w m P mm 

For appropriate matrices P^ v of uj-th roots of unity such that P^\ and P^ are constantly q 
for all li E [m] . 

Proof. For convenience we will assume in the following that 

\Au\ < \Ajj\ for all i < j E [n]. (36) 



This is possible w.l.o.g. by the Per mutability Principle |4.3[ since every matrix A satisfies this 
condition up to symmetric permuta tion s of rows and columns. Let D' = (-D'^) c eZ u with 



L/M = D for all c E Z^. By Lemma 2.1 we have 



%a,d{§, G) = Za,t>'(<P, G) for all digraphs G and all pinnings (f> (37) 



and therefore, EVAL pm (A, D) = EVAL pm (^, D'). By Lemma |4.23| the problem EVAL pm (A, D') 
is #P-hard, if abs (A) has rank at least 2. Assume therefore, that abs (A) has rank 1 that is, 
abs (A) = xy T for non-negative real vectors x and y. The vectors x and y are positive. To 
see this, note that abs (A) is a block, because A is. On the other hand a value Xi = would 
contradict this fact, as then abs (A) i + = Xiy T = implying the decomposability of abs (A). 
Therefore for appropriate w-th roots of unity Qj we have 

Aij = Xiyj ■ Qj for all i,j E [n]. (38) 

By the Hermitianicity of A we have £n E { — 1, 1}. Define q := (u and a diagonal n x n 
matrix n by Ha = q ■ Qi for all i E [n]. Let A' = n^n. The first column of this matrix 
satisfies A' a = HnAnHn = ?CiiCii^j2/iCiiCii = ^iV\ fo r all i £ [n]. As this follows similarly 
for the first row of A' we have 

A' a = q\A' a \ = q\A' u \ = A' u for all % E [n\. (39) 



With D" = (D"[ c l) ceZu defined by D"M = ILD'H the Root Of Unity Transfer Lemma 
yields 



4.7 



ZA,D'{<fi, G) = fn{4>) " -^A',D"(0) G) for all digraphs G and all pinnings 4>. (40) 



5.1 Satisfying Shape Conditions (CI) and (Dl) 
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The function /n is polynomial time computable and non-zero for every </). By definition, all 
diagonal entries of II are u-th roots of unity. Hence IP J_C = II for all c G Z^, that is, D" 
satisfies condition (Dl) . 

We will now perform twin reduction. Let I\, . . . ,I n i be the equivalence classes of the twin 
relation on A' and let r be the twin resolution mapping. By Lemma 4.5 we have 

^4',D"(0> G) = Zyw D [a'j (t o 0, G) for all digraphs G and all pinnings (p. (41) 

Define C = [A'] and D = = (-D [j4 ' 1,Icl ) ce z w . The family D therefore satisfies D\f = 

Dv'v for all i G [n'] and all c G Z w . It follows from the fact that D" satisfies condition 
(Dl) that this is also true for D. Since C is the twin resolvent of A' we have 

Cij = A' for some fj, G ij, v G Ij. 

We may assume w.l.o.g. that 1 G I\ and that for all i,j G [n] and all u. £ £ Ij we have 
| A' | < \A' UU \ iff i < j. This is possible if we assume an appropriate ordering of the I±, . . . , I n i. 
We thus obtain, from equation (39), 

Cn = q\C iX \ = s\C u \ = C u for all i G [n'] (42) 

and, as this was true for A and A', 

\Cu\< \Cjj\ for all i < j G [n'\. (43) 

With rank abs (A) = 1 we have rank abs (C) = 1 and still there are no zero entries in this 
matrix. There are thus positive vectors u and z such that abs (C) = uz T '. By the symmetry of 
abs (C) we have for all i, j G [n'] which implies that u and z are linearly dependent 

with Z{ = — ■ in. Let A = zi • u^ 1 . This rephrases equation (43) to ufX < ujX for all i < j G 

[n'] and therefore the entries of u are ordered increasingly. Let v be the vector of increasingly 
ordered distinct entries of u and w the analogon for z. We have w = \v. Let m be the length 
of w. The definition of C can thus be rephrased as 

/ v lWl P n ... v lWm P lm 
C= : : 

V v mWl P ml . 



for appropriate matrices P^ v of w-th roots of unity. From (42) it follows that P^\ and P-j 1 ^ 



are constant for all \i G [m] — each entry equals Combining equations (37), (40) and (41) 
yields, 

Za,d(<I>, G) = fu{4>) " ^c,d(t" <f>, G) for all digraphs G and all pinnings 0. 
This finishes the proof. n 



Proof (of Lemma 5.5). Assume that EVAL pm (j4, D) is not #P-hard. By Lemma 5.6 we have 
EVAL pin (yl,Z)) = EVAL pin (C, D) for a family D of diagonal matrices satisfying (Dl) and 
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a matrix C which has the following properties. There are vectors v, w G satisfying 
< v\ < . . . < v m and w = Xv for some A > such that 



C = 



P 11 ... viw m P lm \ 



\ v m wiP ml . . . v m w m P mm J 

for appropriate matrices P^ v of w-th roots of unity. By the Hermitianicity of C, there are 
numbers mi, ... , m m such that is an m M x m„ matrix and P' 1 " = P v t l for all E [m]. 
We index C so as to explicate this structure. That is for fi,v G [m], i G [mj and j 6 [m^] 
we let C (/i) j) (l/ j) = \ 
£)I c 1;m ._ £)W 



vPjf ' ■ Index the matrices D' c ' analogously and write, for convenience, 



for all i G [m^]. 



Claim 1. For all h,v,k G [m] and i G [m M ], j G [m„] we have (DM; K , Pj£ • 
ti(D^' K ). Furthermore, equality is assumed only if P?" and P™ are linearly dependent. 



< 



Proof. This is straightforward, we present the calculations for completeness. Application of 
the triangle inequality yields 



fc=i 



fe=l 



tr(I)M ;K ). 



Here, equality can be assumed only if all of the terms Pj k ■ Pf fc are equal 



Claim 2. EVAL pin (C, D) is #P-hard unless the following is true. For all % G [mj and all 
j G [m v ] we have either 

(D^' K ,P™»P£) = for all k G [m] or (D 10 ^,^.^) = tr(D™ ;K ) for all k G [m]. (44) 

In particular the latter case occurs only if P]f* = Pj£- 

Proof. Define, for every q G N the value p = q ■ oo + 1 and a matrix C = C'(p) = C^D^C^. 
For a given digraph G = (V, P), let G p be the result of 2-stretching followed by p-thickening. 
Note that p • d(v) = d(v) (mod oj) for all v G V. Thus we obtain for every pinning <fi 



^c*'(p),d(0 5 G) = Z c ,i)((j),Gp). 



(45) 



The matrix C satisfies 



V V c (p) c (p) p> 101 

Z j Z. j (p,i){K,,k) (v,j)(K,k) (n,k)(K,k) 

K=l fc=l 

m m K 

K=l fe=l 

rrt 



«=1 



5.1 Satisfying Shape Conditions (CI) and (Dl) 
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The definition of p further implies (P^ ) p = Pj* u for all i,j and all p, v. Hence 

m 

= M^ +1 J2^ +1 HD^ K ,P^ . Pp (46) 

K=l 

In particular all diagonal entries of C are positive, as for all p £ [m], i £ [m^] we have 



c. 



.,2(?-w+l) 



^ w 2(^+l) tr(jD [0l; K) _ 



K = l 



If there is a g such that the matrix abs(C') contains a block of rank at least 2, then by 
Lemma 



4.23 



the problem EVAL pm (C, D) is #P-hard. Then further, by equation (45), 
EVAL pm (C,D) will be #P-hard. Assume therefore that, for all q, all blocks of abs(C") 
have rank 1. We shall show that under this assumption C satisfies (CI). 
Consider some 2x2 submatrix of abs (C) of the form 

V \ C (u,j)(p,i)\ \ C (u,j)(uJ)\ J 

By our assumption that for all q this submatrix is not a witness for the existence of a block 
of rank at least 2 in abs (C) we see that either its determinant is zero or C, w a = 0. More 
precisely we arrive at the assumption that, for all p, v G [m], i £ [m(i],j £ [m u ] and 



For all q £ N either C\ 



(ti,i)(v,j) - OT \ C (u,j)(u,j)\\ C (t,,i)M\ 



\ C M(^)\ 2 - ( 47 ) 



If C, 



(fii)(uj) = ^ ^ or infinitely many q then by the definition of C in equation (46), we have 

m 



K=l 



(ni)(vj) = holds only for a finite number of q £ N. Then clearly, 
|Cy jJ 2 must be true for infinitely many such q. By inspection of 



which by Lemma 4.4 implies = (Z)M ;K , PJ£ • P- 1 ^) for all k £ [m] - - as claimed. As- 
sume therefore that ' " 

equation (46), this implies 

m 

£u#«- w+1 >tr(I> 



K=l 



Y,wl {q ' UJ+1) {D^,P™.P l 



i,* I 



K=l 



and thus 



K = l 



K=l 



(48) 



However, application of the triangle inequality yields 



K = l 



< J^wf (D^«,P%.Pp < ^w*tr(uM^). 



K=l 



K=l 
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By C'laimjljwe have 



(£>[°];«,.Pj£ • pp < tv(Dl°^ K ) for each of the k G [to]. Equation ((48) 



therefore implies the existence of some a G U such that 

= wl P (tr(D^ K ) - a(D^ K , P™ . Pp 



K = l 



a(D^ K ,P^.Pp 



As we may choose p arbitrarily large, depending only on £H°J and C, Lemma 4.4 implies that 

tr(D m ' K ) for all k G [m]. 

i „■ ^ — ui ^ iui an rt, in pai u^uiai i ■ ^ — lii 

Lemma 



That is , Pj ^f = olP^I for all ft, in particular Pj\ = olP^ which by our condition from 
that P»\ = P$ implies a = 1. 



5.6 



As C is twin-free, we see that, for every \i G [m] and all i / j £ [m^] there is some re G [m] 
such that Pf* • P/f is not constant. Thus, by the above 

For all /i G [m] and all i ^ j G [m M ] we have (D m '' K , P^ • Pp = for all « G [m]. (49) 

Fix some « and define A K := (D[°1; K )(V2). Let := P^A", then 

= ( D W>", P% • Pp = (T?;, Tp for all i^je [m,]. 

That is, T^ K has full row rank for all /x, k G [m] and by the definition of T^ K it has the same 
rank as P^ K . Thus each P^ K has full row rank and by P^ K = P K ^ it has full column rank, as 
well. Therefore each P^ v is non-singular. Let r = mi be the rank of P 11 therefore P lK being 
a non-singular r x m K matrix for every k G [to] we see that mi = . . . = m m = r. Altogether 

For all fj,, v G [m],P'"' is a non-singular r x r matrix. (50) 

Claim 3. There are mappings ri, . . . r m : [r] — > [r] such that for all n, v G [m] we have 



P.i.i 



P 



for all i,j G [r]. 



Proof. Consider the matrix 



*i 



/ T n ^ 




/ P n A x \ 


, rpml 


1- 


y P^A 1 / 



r 



This matrix has rank r and by the above the first r rows of it form a set of independent 
vectors. In particular, for each \i G [m] and each j G [r] there is an i G [r] such that 

By (T-£,T?l) = (D^\p£»Pll) and Claim|i|we see that Pg = P^ 1 and as C is twin-free we 
further see that for each /i G [m] this is unique. Thus define Tu(i) '■= j and by Hermitianicity 
the claim follows. n 



5.2 The Remaining Conditions (C2), (D2) and (D3) 
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Define a permutation n : [m] x [r] — > [m] x [r] by 7r(/x, i) := (fi, r„ By the above we have 

C( / Lt,i)(i/,i) = v^w u P^ u = v n w vPl\i) Tv {jy Tll us the matrix satisfies 

(C^)( M ,i),(^) = C nMMuJ} = C^ T -i [i))MU)T -i m = VpWvP^-!^^-!^ 

That is 

Define H := ij -1 • P 11 . We have CV^ = ? • Dt« r ® if as desired. 

It remains to prove the technical statements (1) and (2). To see (1), note that (DW^jP^ • 
PP?) translates to (D*, Hj^»H itf ) = H(D^> K )H T . Statement (1) now follows by equation 

We shall show (2), that is, the linear dependence of the Z}M>^. Define the (r — 1) x r 
matrix H' by i^- = HijH r j for all j G [r] and i G [r — 1] and let diag(£)M'^) be the vector 

of diagonal entries of then (1) implies that (HD^'^H T ) r j = for all j ^ r and all 

G [to]. This is equivalent to H' diag(D^'^) = for all fj, € [to]. As H is non-singular we 
have rank H' = r — 1 and therefore all solutions to this system of linear equations are linearly 
dependent. This proves (2). n 

5.2 The Remaining Conditions (C2), (D2) and (D3) 

We now turn to the part of the proof devoted to the remaining conditions (C2), (D2) and 
(D3). 

Lemma 5.7. Let C, D satisfy conditions (CI) and (Dl) then either EVAL P '"(C, D) is #P- 
hard or conditions (C2) and (D2) hold. 

Lemma 5.8. Let C, D sate/y conditions (C1),(C2) and (Dl) and (D2) then either EVAL P "'(C, D) 
is #P-hard or condition (D3) holds. 



The proofs of both lemmas are technically similar. Therefore we start by discussing the 
technical details they have in common. The most technical part is given by the Pre-Uniform 
Diagonal Lemma |5.9| below. Further, as the #P-hardness proofs in both, Lemmas 5.7 and 



5.8 are based on the same reduction, we discuss this reduction in Lemma 5.10 



A Technical Tool. We call a diagonal-matrix d G C mxm pre-uniform if there is a d G C 
such that for all i G [to] we have Da G {0, d}. 

Lemma 5.9 (Pre-Uniform Diagonal Lemma). Let H G U£, xr be a non-singular matrix 
and D G IK^ ><r a non-negative diagonal matrix. Let K = {k G [r] | ^ 0} be the set of 
non-zero diagonal entries in D then one of the following is true 

• There is a computable p$ G N such that for allp > po the matrix abs (^HD p H T ^j contains 
a block of rank at least 2. 

• D is pre-uniform and for all I C [r] the following holds. If Hjk is non-singular then it 
is a complex Hadamard matrix. 
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Proof. The matrix B := HD P H satisfies, for all i,j G [r]: 



B. 



^2 H ikH jk D P kk 



k=l 



^ HikHjk 
k&K 



D p 
u kk 



That is, for every / C [r] we have Bn = HikD p kk H IK . Let k = \K\ and fix a set / C [r] 
such that |/| = k and -f/j/f has rank k. Note that every principal 2x2 submatrix of abs (Bji) 
has non-zero determinant. To see this, fix an arbitrary such submatrix induced by distinct 
indices i,j G [r]. Define a matrix X = Hj^D^^ which is non-singular since Hjk is. We 

— T 

have B = XX . The determinant of the submatrix just defined is thus 

\(Xi : *, Xi^)\ | (Xi^, Xj^} 



(H* )K D p KK H^ K )ij. 



(51) 



\(Xj^,Xi^)\ | (Xi t *, Xj^}\ 



\(Xi^, Xi^)\\(Xj^, Xj^) \ — I (Xi^, Xj^) Y 



Using the Cauchy-Schwarz inequality, the non-singularity of X implies that the above deter- 
minant is non-zero since 

Therefore, every principal 2x2 submatrix of abs (Bjj) has non-zero determinant. The ex- 
istence of a block of rank at least two in abs (Bjj) is thus guaranteed, if we find a non-zero 
off-diagonal entry. Hence, to finish the proof it remains to show the following: If there is no 
computable po such that for all p > po the matrix Bn contains a non-zero off-diagonal entry, 
then D is pre-uniform and Hjk complex Hadamard. 

For every pair of values i,j G / define a family Z^ of equivalence classes Z G Z^ such 
that k,k' G Z if, and only if H^Hjk = Hi^Hj^i. Define for each such Z G Zij the value 
£z '■= HikHjk for some k G Z. Furthermore, define a family J of equivalence classes induced 
by the equality of diagonal entries of Dkk- That is, for each J G J we have j,j' G J iff 
and denote by dj the corresponding diagonal value of D. We obtain, for all 



D j'f 



keK 



ofj tz\Jnz\ 



By Lemma 4.4 there is some po such that the following holds. For all p > po if Bij is zero 



then, for every J € J, 







y] (z\j^z\ = y~] H ik Hj k 



zeZi 



fceJ 



As this is true for alH 7^ j G I the matrix Hjj is non-singular and complex Hadamard. This 
implies that J = K, that is, D is pre-uniform. n 



The #P-hardness construction. 

Lemma 5.10. Let C G C™ xm and D a family of diagonal matrices which satisfy (CI) and 
(Dl). For allp,q G N, there is a diagonal matrix 



J2v p D^' K 



K = l 



2q 



(52) 



such that the following is true. 

If abs (HQH I contains a block of rank at least 2 then EVAL P '"(C, D) is j^V-hard. 



5.2 The Remaining Conditions (C2), (D2) and (D3) 
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Figure 1: The reduction template T P)q for p = 3, q = 1. 



Proof. We define a reduction template T PtQ which will be used in the following reductions. 
This template is a digraph T Pt q = (T^g, E Pt q) with two terminal vertices u and t> which are 
connected by a length 2 path with middle vertex w. The most important subgraph of this 
template is a graph P(p), which consists of two vertices b and a connected by p many paths 
of length 2 directed from b to a. Then, q many disjoint copies of P(p) are attached to w 
by identifying their terminal vertex a with w and further q copies of P(p) are attached by 
identifying their terminal vertex b with w. Figure [T] illustrates the construction. Its formal 
definition is given by 

Vp, q = {u,v,w,Xi,Xij,yij,yi \ i G [q],j G [p]} 

E Ptq = {uw,wv,XiXij,XijW,wyij,yijyi \ i G [q],j G \p]} 

For a given digraph G, let G' be the graph obtained from G by replacing every edge uv by a 
distinct copy of T Piq . Let 4> be a pinning of G, then 

^ jD (0,G) = 2 c ,d(^G !/ ) (53) 

where C' is a matrix whose structure we will analyze in the following. We start the analysis 
by observing a technical detail. We have, by Lemma |5.5| (1), 

m m 

CD^C = (vv T ) ®^2wlHD^H T = (vv T ) <g> ^v^\x{p^)I r 

v=\ u=l 

Condition (Dl) implies that tr(D^' u ) is positive for all v and thus there is some positive 7 
such that 

C£>[°JC = (vv T )®jI r . (54) 
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Here I r denotes the r x r identity matrix. Denote by Zcj)(i,j;T pq ) the partition function 
of T p>q with vertex u pinned to i and v pinned to j. We have C'^ = Zc,T){i,j]T P;q ) for all 
Further the definition of T Pi g directly implies that C' = CAC for some diagonal matrix 
A = A(p,q). In the following we will index A by pairs (/U,i) £ [to] x [r] and further define 



A ii = a (aM)(aM)- We obtain 



\k=1 k=l 



1 / m r 

I EE4j I;K H ^ 

\k=1 k=l 



(K,k),(fi,i) 



Applying equation (54) we see that, for all (n, k), £ [to] x [r] 



CD m c 



(/x,i),(/c,fc) 



CD n c 



(K,k),(fi,i) 



which is zero unless k = i. By condition (Dl) we have Z)I-p1 = DM and thus we can rewrite 



r 



\K=1 



\K = 1 



1 / m 



e^? k 



\K = 1 



VK=1 



And thus 



K=l 



2q 



A direct computation yields 



(55) 



CAC rr J 8)iJ | E^^ | H ' 



(56) 



and by equation ( 55 ) we obtain 



fi=i fi=i 



K = l 



2q 



■l 2pq ■ G 



where the last equality follows from the definition of = Q(p, q) in equation (52). Equation 



(56) then rephrases to 



CAC = vv T ®H ( 7 2p,? G) H T = 7 2w • vv T <g> HQH T . 



It thus follows that abs (CAC) has a block of rank at least 2 iff abs yHQH j does. Assume 

that ab s ^H@ H T ^ contains a block of rank at least 2. Then EVAL pin (C / , D) is #P-hard by 

Lemma [423| Equation ([53]) implies that EVAL pin (C, D) < EVAL pin (C, D). This finishes the 
proof. n 



5.2 The Remaining Conditions (C2), (D2) and (D3) 
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The proofs of Lemma 5.7 and 5.8[ 

Proof (of Lemma\5/fy. Assume that EVAL pin (C, D) is not #P-hard, we shall show that (C2) 
and (D2) are satisfied. Fix a posi tive p such that [p] = 0. Let = 0(p, q) the matrix defined 



5.10 



for all p,g£N the matrix abs (HQH T ) contains no blocks 



in equation (52). By Lemma 
of rank at least 2. Recall the definition of = Q(p, q) and note that our choice of p yields, 

|2g \ 2q 



K = l 



E«^ I01; iE<^ I01;K 



\K = 1 



The second equal ity follows by condition (Dl) as the absolute values | • | are inessential for 
DM. By Lemma 



5.5 2), we have 



2q 



\K=1 



Define f(p, q) := Y?u=i w fi v fi Pq dfi {Y^=i v K.d K ) zq which is positive for all p,g£N. Then 

m / m \ 2g 

©= E^X P % E-x (^ I0I;1 ) 2g+1 = /(^)P I01;1 ) 29+1 - 

^t=l Vk=1 / 

For notational convenience let D := DM; 1 . Then 



HQH T = H(f(p, q) ■ D 2q+l )H T = f(p, q)HD 2q+1 H 



which implies that abs ^HQH T J contains a block of rank at least 2 iff abs ^HD 2q+1 H 
does. 

We therefore arrive at the assumption that for all q G N the matrix abs (^HD 2q+1 H 
does not contain a block of row rank at least 2. The Pre-Uniform Diagonal Lemma |5.9| thus 
implies that H is complex Hadamard and D = al r for some a. This proves (C2). Let 
AjiJJ = d^a with the as given by Lemma 5.5 '2). Then L>M = A^' I r which proves 
(D2). □ 



tT 



Proof (of Lemma JM. Assume that EVAL pm (C, D) is not #P-hard, we shall show that (D3) 
is satisfied for for all c G Z^. Fix an arbitrary c G Z w and let T c = {p G N | p = 
c (mod cj)} be the set of natural numbers which are congruent to c modulo oj. That is, we 
have \p\ = c for all p G T c . Let = 0(p, g) be the matrix defined in equation (52). By 
Lemma 15.101 we have 

For all p, q G N the matrix abs [HQH] contains no blocks of rank at least 2. 



We shall see how this implies condition (D3) for c. Assume from now on that p G r c . Then 
by the definition of = Q(p, q) in equation ( J52| ) we have 

n 2 ? 



11=1 



E^ bI; 

K = l 

n 



K = l 



2q 
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where the last equality follows from [pj = c = [cj and our assumption that condition (D2) 
is satisfied. Let f{x) := (EJ£=i w^aJJJ) and define 



A 



K=l 



(57) 



Then 6 = /(2pg) • A 2 <? and we have #6F" = H(f(2pq) ■ A 2c i)H T = f(2pq)HA 2q H ± . That 
is, the existence of a block of rank 2 in abs (^HQH T ^j is equivalent to the existence of such a 
block in abs (HA 2q H T ^) and thus it follows from our assumption that 

For all p,g£N the matrix abs ( H A 2q H T ) does not contain a block of rank at least 2. 



(58) 

Claim 1. For all p G T c the matrix A = A(p) is pre-uniform. 

Proof. Assume otherwise. Th at is , fix a p G r c such that A(p) is not pre-uniform. Then the 

implies that there is a q such that abs (^HA 2q H T ^j contains 

H 



5.9 



Pre-Uniform Diagonal Lemma 
a block of rank at least 2, in contradiction to equation (|58). 



In the above claim, we used the parameter q to show that A is pre-uniform irrespective of 
our choice of p G T c . In the following we will show that choosing p large enough we can infer 
condition (D3) for c. To achieve this, we will prove the following two statements 



Claim 2. There is a po such that for all p G T c with p > po, 

If Ajj = A(p)a = then for all \i G [m] we have Djf"^ 



0. 



Proof. Fix an i G [r]. By the definition of A in equation (57) we see that An = if, and only 
if, 



K = l 



If this equation is satisfied for all p, Lemma 4.4 implies the existence of values pi such that 
for all p G r c with p > pi we have jD|? = for all fi G [m]. We derive such a p» for each 
i G [r] and the claim follows for po = maxjpi, . . . ,p r }. H 

Claim 3. There is & p= G N such that for all p £ T c with p > p = and all i,j £ [r] we have 
that 

An = Ajj implies that there is a Qj G U such that = CijDj^'^ for all \i G [m]. 

Proof. By the definition of A the assumption An = Ajj is equivalent to 



[c];k 



[c];n 



K = l 



K=l 



5.3 Finishing the Proof of Lemma 



5.1 
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Or likewise, there is a Cij G U such that 



m 

K = l 



By Lemma 4.4 there is a such that for all p > p™, the above implies that all coefficients 
d|j C "' k — djD^j' K are zero. The claim follows with p= = max{pij \ i ^ j G [r]}. _| 

With these three claims we are now able to finish the proof of the lemma. Let p G T c such 
that p > max{po,p=} with po and p= as in Claims [2] and |3j We will argue by considering the 
diagonal entries of A = A(p). Define M = {// G [m] | / 0} and let A = {i G [r] | A« / 

o}. 

If A = then Claim [2] implies that =0 and the proof follows with A' c ' =0 and 
XjM = 0. 

Assume therefore that A / and fix some a G A. As A is pre-uniform by Claim [I] we 
have Ajj = A aa for all i 6 A. By Claim [3] this implies that 

Dlf' M = Cw-DMf for all M G AT, * G A. 

Define A@ = Z# M for every ^ and note that Z?|a ^ = if, and only if p, £ M. Define 
U\f = Qa for all i G [r]. Then ijW = AW (8) f/W. □ 



5.3 Finishing the Proof of Lemma 5.1 



Let A be a connected non-bipartite w-algebraic Hermitian matrix and D a diagonal matrix 
of positive vertex weights. By Lemma 5.5 either EVAL pm (A, D) is #P-hard or the following 



holds. There is a matrix C satisfying condition (CI) and a family D of diagonal matrices 
which satisfies condition (Dl) such that, 



Further, by Lemmas 



5.7 



and 



EVAL pin (A, D) = EVAL pin (C, D). 

the problem EVAL pm (C, D) either is #P-hard or the matrix 



5.8 



C satisfies conditions (C1)-(C2) and the family D satisfies (D1)-(D3). By Lemma 5.3 
then have 

EVAL pin (C,D) = EVAL pin (#,U) 
where H and U define an (H— STD) problem. This finishes the proof. 



we 



6 Hadamard Components 

In this section we will prove Theorem |3.5| by proving the following two lemmas. 

Lemma 6.1. Let H be an u-algebraic n x n matrix and D a family of diagonal matrices 
defining an (H STD) -problem. Then either EVAL p '"(ff, D) is #P-hard or H andT) satisfy 
conditions (GC) ; (R1) through (R5) and the Affinity Condition (AF). 

Lemma 6.2. Let A be u-algebraic with underlying nxn block H and D a family of diagonal 
matrices defining a (B H— STT>) -problem. Then either EVAL P '"(A, D) is #P-hard or H and 
D satisfy conditions (GC),(B Rl) through (B R5) and the Affinity Condition (B— AF). 
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Figure 2: The gadget for q = 1 



Proof (of Theorem 3.5). Let A be an u;-algebraic matrix and D a family of diagonal matrices 
denning either an (H— STD) problem or a (B H STD) problem. 

If EVAL pin (A, D) is not #P-hard then the following holds. If A and D define an (H STD) 



problem then Lemma 6.1 implies that A and D satisfy conditions (GC), (Rl) through (R5) 
and the Affinity Condition (AF). If A and D define a (B H STD) problem then Lemma 6.2 
implies that A and D satisfy conditions (GC), (B Rl) through (B R5) and the Affinity 
Condition (B— AF). This finishes the proof of the Theorem. n 



Throughout the section, n and u will have a fixed meaning. That is H will be an w-algebraic 
nxn matrix. We will discuss the non-bipartite case in order to prove Lemma 6.1 In a sequence 
of several steps the proof will show that the problem EVAL pm (ff, D) is #P-hard unless the 
conditions (GC), (Rl) through (R5) and the Affinity Condition (AF) are satisfied. This 
will be done in an inductive manner always relying on the conditions for which we have 
shown this so far. The proof of Lemma |6.2| follows by similar means, therefore we omit it. 
The interested reader may find the proof in |Thu09j . 



6.1 The Group Condition (GC) 

The first step of the proof is to enable a group theoretic description of the matrix H underlying 
the (H STD) problem at hand. 

Lemma 6.3. Let H be a matrix and D a family of diagonal matrices defining an (H STD)- 
problem. If H does not satisfy the group condition (GC) then ~EVAL pm (H, D) is f^P-hard. 

Proof. Let G = (V, E) be a given digraph and 4> a pinning. Let G q = (V^, E q ) be the digraph 
obtained from G by 

V q := V U {v e , v' e , v e>1 Ve t2q \ e £ E} 

E q := { uv e> 2i, v e ,2i-iu, v e> 2iV, vv e ,2i-i | i € [q], e = uv £ E} 

U{v e V e ,2i, V e ,2i-lV e , V e ,2iv' e , v' e V e ,2i-l \ i € [<?], e £ E} 

We can think of G q being obtained from G by replacing every edge by a distinct copy of a 
graph T q which, for q = 1 is illustrated in Figures] Let i,j, a, b be spins assigned to u, v, v e ,v' e . 



6.1 The Group Condition (GC) 
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This witnesses a reduction EVAL pin (CM) < EVAL pin (F, D) with 

n n / n \ 1 / n \ 1 

c]j' = E E ( E ^ ic ^ c ^ ac ^ cb j I ^ H c jHj c H ca Hbc J 

a=l 6=1 \c=l / \c=l / 

n n 

= e e K-^a,* * Hf,^,Hj^ • Hj^)\ 2q 

a=l 6=1 

By an argument based on c]| for i,j G [n] we will see in the following that EVAL pm (C^) is 
#P-hard if H does not satisfy the group condition. We will not give a direct proof, but we 
will rather show that, 

if there are a', b' G [n] such that H a >^ • Hy^ then EVAL pin (#, D) is #P-hard. 

To see that this finishes our proof, we shall prove that 
Claim 1. If Hi^ • Hj } * G R{H) for all i, j G [n] then H satisfies (GC). 

Proof. If Hi t * • IZ^* G R{H) for all t, j, G [n], then there is a mapping r : [n] 2 — > [n] such 
that, for all i,j we have Hi * • ii^* = H^jy^. Fix an arbitrary j G [n] and consider the map 
Tj := r(-,j). This map is bijective. To see this, assume that r(a,j) = r(c, j), which implies 
that H a ^ • Hj^ = H c ^ • i?^* which yields H a = H c . Therefore, we see that, for all i,j G [n] 
there is a c G [n] such that r(c,j) = i, that is, H c ^»Hj^ = H^* which yields Hi^oHj^ = H c ^. 
Since H is Hermitian this also proves the second part of the group condition. -| 

Assume that there are a' , b' G [n] such that H a /^ • Hy* ^ R(H). Fix such a' , b' and let i = 1 
then we have 



c'fj — ^2 E K^ a >* • Hb,*iHj,* • H\,*)\ 2q — E E l^ a '* * Hb,*, Hj^)\ 2q 

a=l 6=1 a=l 6=1 

Note that for all a, b we have \{H a ^ • -f^,*, ^j,*)! < n and further 
Claim 2. \{H a i^ • Hy^,Hj^}\ < n for all j G [n]. 

Proof. Assume that \(H a /^ • Hy *, = re. Then the Cauchy-Schwarz inequality implies 

that there is a ( such that i? a ',* • Hy ^ = C-^j,* but as H is normalized we have Hji = 
(H a /^ • Hy^)\ = 1 and therefore £ = 1 which implies -ff a ',* • Hy * G R(H) in contradiction to 
our assumption. -I 

Claim 3. There is a j G [n] such that {H a /^ • Hy Hj *) ^ 0. 

Proof. Otherwise, H a /^»Hy^ would be orthogonal to all vectors in R(H) and by H a i *»Hy ^ 
R{H) this would imply an n + 1 element basis of an n-dimensional vector space. h 

By construction we have abs (C™) = C™ and is symmetric for all q. Therefore the proof 
follows from Lemma 4.23 if we can show that there is a j G [n] such that the 2x2 submatrix 
of induced by 1, j is indecomposable and has rank 2. We have c\f = n 2q+l for all i G [n], 
as 

n n n n 

cf = E E K ff ^ • ^,*^,* • ^,*>i 29 = E Ek^.*>^*>i 2<2 

a=l 6=1 a=l 6=1 
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and the term \(H a ^,Hb^)\ is non-zero iff a = b. Fix j G [n] as in Claim [3j We see that the 
2x2 submatrix of induced by 1, j satisfies 

/ Cg eg \ = / n 2 ^ 1 eg \ 

I ^ ^ J " I ) ■ 

Claim [3] further implies that it is indecomposable. To show that this submatrix is a witness 
of a block of rank at least 2 it remains to show that its determinant n 4q+2 — (Cg) 2 is non-zero 
for some q. Define values Co, . . .c n by c v = |{(a, 6) G [re] 2 : \(H a ,* • Hb^,Hj^}\ = Then 
we can rewrite 

n n n 



eg = X) i ( Ha '* • Hb >*> H i,*)\ 2q = Y, c ^ 2q 

a=l b=l u=l 

Assume, for contradiction that the above mentioned determinant is zero for all q. This implies 
that 



u=0 



for c' , . . . , c' n such that c' = c' n _ 1 = and c' n = n. By Lemma 4.4 we see that for some large 
enough computable q this implies that c v = c' u for all v E [n\. In particular n = c' n = c n 
which contradicts the fact that Claim [2] implies c n = 0. n 

6.2 The Representation Conditions (Rl) through (R5) 

Now we shall describe how to satisfy the Representation Conditions (R1)-(R5). 

Lemma 6.4. Let EVAL p/ "(.H~, D) be an (H-STD) problem such that H and~D satisfy (Rl)- 
(R2). Then ~EVAL pm (H, D) is #P-hard unless conditions (R3) and (R4) are satisfied. 

Proof. Fix some c E X u , we may assume that c > as for c = we have /3q = and <5q = & 
by the fact that L>M = j n . Further, for -c> we know that D M = by (H-STD) 

and thus A c = A_ c . 

For any c > we will construct a reduction as follows. Let G = (V, E) be a given digraph 
and (j) a pinning. We construct a digraph G c from G by replacing each edge e = uv by a 
distinct length 2 path and call the middle vertex of this path v e . Then we add another vertex 
v' e to the graph and connect v e to v' e (in this direction) by c distinct length 2 paths. Formally, 
we have G c = (V c , E c ) with 

V c = V U {v e ,v' e ,v etl , . . . ,v etC \ e G E} 
E c = {v e v e ^,v et iv' e I i G [c], e£B) 

We have c = <9(t> e ) = — 3(i>e) and therefore the vertex weights of v e and t> e have the same 
support A c . The above construction witnesses a reduction EVAL pin (C, D) < EVAL pin ( J ff, D) 
for a matrix C which we will describe now. We label the rows and columns of C symmetrically 
by the elements of (5. Let a,b, g, h G (5 denote the spins of u, v, v e and v' e respectively. Then 



c 



C a , b = Y, Ba,- 9 H 9 ,- b DS DjTfl J2 H^H^Dl^ 
gee he® V/^ee 



g,h£A c 



6.2 The Representation Conditions (Rl) through (R5) 
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The second equality follows from the fact that = I n and H^—h = Hh—u,- We know that 
{H g) », -H/i,*) = unless g = h. And if g = h, then (H g ^, H^) = n. Thus, we have 



c a , b = E 

ff6A c 



n 



9 eA c 



Note that the right equality uses H a - g = H a>g . We know that A c C & but it might be the 
case that the neutral element of is not contained in A c . Fix some A € A c and define 
S = {g — A | g G A c } which clearly contains 0. We have 



Cab = n C Hb-a,g = n° Hb- a ,x+\ = n C Hb- a ,\ ■ E Hb- a ,x 

geA c x&S xeS 

Let x : © - ► {0, 1} denote the characteristic function of 5, then row of C satisfies 

C ob = n c H b ,\ E H b>x x(x) = n c H b: x(H b ^,x}- 
xe® 

Define C = {b G © | (6,fli) w = for all 5 G 5} and recall that [0] = 1. Thus 

<H 6 ,*, X ) = E X ^ -9)-] = E ^ HMM =£«-[0] = |5| for all 6 G C . 

9 e5 9 eS 9 es 

Claim 1 If EVAL(C, D) is not #P-hard then for all b G & \ C we have (-H& *,x) = 0. 



(59) 



(60) 



Proof. Fix an element b G (5 \ Cq. Equation (59) straightforwardly implies that C aa = n c \S\ 



for all a G 0. Consider a 2 x 2 submatrix of abs (C) induced by and 6, 



I Coo | \Cob\ 
\Cw\ \C bb \ 



n 



\S\ \(H b ,„ X )\ 
\(H b ,*,x)\ \S\ 



By Lemma 



4.23 



the problem EVAL pm (C, D) is #P-hard, if abs (C) contains a block of row 
rank at least 2. Therefore no such 2x2 submatrix can be a witness for the existence of a 
block of rank at least 2 in abs(C). This however is possible only if either \{H b ^,x)\ = or 
the determinant of this submatrix is zero. The latter would imply \S\ = |(.H&*)X)I which 
is impossible for the following reasons. First of all, as H is normalized, and G S the 
expression (-E/&,*, x) = J2xes ^b,x contains at least one 1 entry. However, as b £ Cq we see 
that Hb^ = 3£ [(b, *) u ] is not constantly 1 on S which implies \ (H b ^,x)\ < \S\ and the claim 

follows. -| 

Define 6 as the smallest subgroup of & containing S. Recall that by the subgroup criterion 
of finite groups, each element of & is of the form S 5 es A 9 g for appropriate A 9 G Z. Thus for 

every b G Cq and every g' G 6 we have (b,g') u = {b, (j2 g eS ^g9j)u = E g6 s A <?(&> 9)u = °- 

Let x& '■ ® ~~ ^ {0> 1} De the characteristic function of 6. Analogously to the above we see 
that for all b G Cq, 



(H b ,*,xe)= [-<&,</%]= 

g'ee g'ee 



EVM 

ges 



E^[0]= |6|- (61) 
g'ee 



We will furthermore see that 
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Claim 2. For all b G & \ Co we have (-H&,*, xs) = 0. 

Before we prove this claim we will argue that this finishes the proof of the lemma. Note first 
that, the basic properties of Fourier analysis of Abelian groups imply that 



and 



X&{x) 



jgjE^*' 



Xe) 



H 



9,x- 



9 e« 



Combining Claims [T] and [2] and equations (60) and (61) we see that for all g G (5 we have 
{Hg t *,xe)\S\ = (H g ,*, x)|©|- Which implies that = Xe^l^l for all x G &, that is 

(3 = 5. It follows straightforwardly that A c = A + 6. That is, it is a coset of 6. This 
proves condition (R3). Recall then from the discussion of conditions (Rl) through (R5) in 



Section 3.3 that condition (R4) is satisfied, as well. This finishes the proof. 



The Proof of Claim [2} Let b G 6 \ C and define G b a = {g G 6 | (b, g) w = a} for all 

a G As (&, *) u defines a homomorphism, each nonempty set & b a is just a fiber of this 
homomorphism and Sq is its kernel. In particular, Sq is a subgroup of 6, it is non-empty as 
G ©o an d the other non-empty sets G c a are its cosets. 

Let 21 = {a G Z w | G b a 7^ 0} be the subset of Z w which selects all non-empty 6^. The set 
21 is an Abelian group. To see this, recall that G 21 and further the existence of elements 
g G 6 b a and h G & b a , denotes (b, g) u = a and (b, h) u = a'. We have g + h G & and by the 



bilinearity of the (=1 



operator a + a' = (b,g) u + (b, h) w = (b,g + h) L 



non-empty and therefore the subgroup criterion implies that 21 is a group. We have 



Thus ® b a+a , is 



{Hb,*,Xe) 



E 

gee 



ar[(b-g) 



E E 



6° 



«e2t #e6 

(fe,9>uj=-o 



ae2l 



As the (3q are cosets of 6q they all have the same cardinality and we obtain 



(Hb,*,Xe) 



axbl 



E^m 



It remains to show that the right hand side sum is zero. By the Fundamental Theorem of 
Finitely Generated Abelian Groups, there is a direct sum decomposition 21 = <t\ @ . . . @ <t z 
and each of these cyclic groups has all elements of the form Aj/ij for some Aj G ^ id(hi) an d 
hi for i £ [z]. Therefore, 



ae2t 



a 



, E 

Aiez ord ( hl ) 



E 



i=l 



And each of the sums on the right hand side satisfies 

E ^N A * : 



E 

8—1 Aiez ord ( hi ) 



(62) 



6.3 The Affinity Condition (AF) 
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provided that hi ^ 0, as [hi] is some ord(^)-th root of unity. Further, since b G <5\Co there 
is at least one element in 21 which is not the neutral element. Therefore in the decomposition 
of 21 at least one non-trivial generator hi ^ exists. This proves that the term in equation 



( 62 ) is zero and thus finishes the proof of the claim. n 



Lemma 6.5. Let EVAL pin (H, D) be an (H STD) problem such that H satisfies (GC). 
Assume that H and D have a representation as given in (R1)-(R4). There is an (H 
STD) -problem KVAL pin (H, D') which satisfies conditions (R1)-(R5) such that 

EVAL pin (H, D) = EVAL^^D'). 

Proof. Define, for each c G Z w values A c as follows. If = let A c = 0. Otherwise, 
we know by condition (R4) that = JT"[p c (0)]. We define A c = [/^(O)]- 1 and let 

D' = (A c • ZH^) c gz w . Let G = (V, £7) be a digraph and ^ a pinning. Let, for each c G Z w be 
n c the number of vertices v £ V \ def(4>) such that [9(u)J = c. Then 

This proves the reducibility in both directions. n 



6.3 The Affinity Condition (AF) 

In the present situation we are faced with an (H— STD) problem EVAL pm (iJ, D) which 
already satisfies the group condition (GC) and the represent ability conditions (R1)-(R5). 
Our aim will be now to provide the last big step necessary for proving Lemma |6.1| This step 
is given by the following Lemma. 

Lemma 6.6 (The Non-Bipartite Affinity Lemma). Let EVAL P '"(F,D) be some (H- 
STD) -problem such that H satisfies (GC). Assume that H and D have a representation as 
given in (R1)-(R5). // (AF) is not satisfied then EVAL p '"(i/, D) is #P-hard. 



The proof relies on a construction which we will analyze separately in Lemma 6.7 so as to 
make it more digestible. We shall describe the reduction first. 



The Nymphaea Reduction. We consider the problem EVAL pm (H , D) as above. Let 
G = (V, E) be a given digraph and <j) a pinning. For parameters p and q, we construct a 
digraph G' = G'(p, q) by replacing each edge in G by a graph T Ptq which has two distinguished 
terminals u and v. We refer to this graph as the "Nymphaea" (see Figure [3]). The vertex set 
of T m is 

{u,v,z,w,Uij,Uij,Ui,Ui,Vij,Vij,Vi,Vi,Xi,Xi,yi,yi \ i £ [p],j G [q]} (63) 
Its edge set is 

{uuij, ViUij, UijUi, UijXi | i G [p], j G [q]}U {uiju, UijVi, umj, XiUij \ i G [p], j G [q]} 

U {vijv, Vijm, ViVij, xmj | i G [p], j G [q]}U {vvi,j, ym,j, v~i,jVi, v~i,jXi \ i G [p], j G [q]} 

U {xiz, wxi, zyi, yiW \ i G [p], j G [q]}U {zxi, XiW, y~iZ, wyi \ i G [p], j G [q]} 

(64) 
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Figure 3: Nymphaea for p = 1, q = 1 



Lemma 6.7. Let EVAL pin (H, D) be an (H-STD) problem such that H satisfies (GC). 
Assume that H and D have a representation as given in (R1)-(R4). 
Then for all p G N and q G Z w £/te Nymphaea reduction witnesses 

EVAir n (C) < EVAir n (F,D) 

/or a non-negative real valued n x n matrix C = C(p, q) which satisfies the following. For all 
u, v £ f3 g + & q we have 



Cu,v — ^ 

gee 



^2 ^ M V - Pq + X )~ Pq( U ~ Pq + X ) + (ffi X ) 



2/J 



Proof Let G = (V, i?) be a digraph and (/> a pinning. For parameters p and q let G" = G'(p, q) 
by the digraph of the Nymphaea reduction as described above. 

Observe first that there is a certain axis of symmetry in the Nymphaea which has also been 
made explicit in Figure [3j The vertices u,v,z,w lie on this axis and all other vertices occur 
in two forms, without the bar ( e.g. Ui) and with the bar {ui). Further, the connections of 
the unbared vertices are the same as for the bared ones except for the fact that the directions 
are reversed. By H being Hermitian, this amounts to a conjugation of the corresponding 
contributions and will become important in the following. First of all, due to this symmetry, 
it will suffice to analyze the unbared part of the gadget. The results for the bared part then 
follow clearly by conjugation. 



6.3 The Affinity Condition (AF) 
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Analyzing the Nymphaea. To simplify the analysis further we will focus on one of its sub- 
graphs first. For some i G [p], j G [q] consider the subgraph induced by vertices u, Uij, Ui,Xi,yi. 
We will fix the spins of u,Ui,Xi, and m for the moment and — slightly abusing notation - 
denote these by just the vertices they correspond to. We denote the contribution of this sub- 
graph excluding the vertex weights of u, Ui,Xi, yi by Z\{u, Ui,Xi, yi). The vertex mj will not 
be fixed and its spin will be denoted by ji. Note that d(uij) = which implies by (H— STD) 
that the vertex weights on Uij are inessential. Then — recall that H is indexed by elements 
of the group © — it follows that 

Zi(u, Ui, Xi,yi) = H Ut —f t Hf it — Xi Hy it —f l Hf it — Ui 

= ^2 H u - tl H Xi - fl H yi - tl H Ui - fl 
MS© 

This can be rephrased as Zi(u,m,Xi,yi) = (H u - Xi+yi - Uii *, Ho,*) which by the Hadamard 
property of H implies that Z\(u, Ui, Xi,yi) = unless u — x% + yi — U{ = and if so, we have 
Zi(u,Ui,Xi,yi) = n. 

Analogous reasoning on the subgraph induced by vertices v,Vij,Vi,Xi,yi shows that its 
contribution is 

Z2(v,Vi,Xi,yi) = (H Vi - v+Xi - yu *,Ho t *). 

Similarly, this is zero unless v j — v + Xj — yi = 0, and if so we have Z 2 (v, Vi, Xi, yi) = n. 

Let now Z q (u, Ui,Xi, yi) denote the contribution of the subgraph defined by u, Uij,Ui, Xj, yi 
for all j G [q]. Then straightforwardly Zf(u, Ui, Xj, y{) = Zi(u,Ui,Xi,yi) q . Defining Z q (v, Vi, Xj, yi) 
analogously on v, v itj ,Vi, x i: yi for all j G [q] we obtain Z%(v,Vi,Xi,yi) = Z 2 (v,Vi, Xi,yi) q . Fur- 
thermore, for the subgraph on vertices z,w,Xi,yi we have a contribution of 

H X i,—zH Zt —y i H Wt —x i Hy i — w = H Xi — z Hy i ^ z H Xi — w Hy i — w = 3s [(xj — yi,z w) tJ ] 

Let us now turn to the complete gadget. Note that = d(z) = d(w) = and <9(xj) = d{yi) = 
for all i G [q\. Therefore the corresponding vertex weights are inessential. The remaining 
vertices Ui, Vi satisfy —d(ui) = d(vi) = q. By condition (H STD) we know that = DM 

which particularly implies that the support of these two is the same, namely A q = A- q . Let 
fj,, v be the spins of Ui and Vi respectively. For every fixed i G \p] and every configuration 
of u,v,Xi,yi the subgraph defined by u,v,Xi,yi,Ui,Vi,Ui t j,Vij for all j G [q] then yields a 
contribution 



Z(u,v,Xi,yi) = ^2 D^j l Dl q lZ 1 (u,^,Xi,yi) q Z 2 (v,iy,Xi,yi) q 

and by the above reasoning on Z\ and Z 2 we have 

Z(u,v, Xi , yi ) = YI 4!i^l(^-^+«-M,*^o,*) 9 (^-, +:Ci - 2/i> *,^o,*) 9 (65) 

= n 2q D M D M (66) 

11 ^u-Xi+yijU-Xi+yi-^v-Xi+yijV-Xi+yi \ UU J 
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where we now included all vertex weights. By the fact that the bared part of the Nymphaea 
has all edges in reverse direction we obtain by the same arguments as above that 



Hz ! —xiHy it —zH Xi! — w H Wi —y i — X \(xi Hi, z w)^] 

Further the the subgraph defined by u,v, Xi,y~i,Ui,Vi,Uij,Vij for all j £ [q] then yields a 
contribution 



Z(u,v,Xi,yi) 



That is 



Z(u, v, Xi,yi) = Z(u, v, Xi, yi) (67) 

Note further, that by construction all original vertices have grade 0. Altogether, the Nymphaea 
therefore witnesses a reduction EVAL pm (C / ) < ~EVAL pm (H, D) for a symmetric non- negative 
matrix C defined by 

c 'u,v = III X! Z{u,v,Xi,yi)X [{xi - y^z - 

z,w£<& i=l \xi,yi£<8 



X Y{\ X z {u,v,Xi,yi)X [(xi 

Hit Z W/^l 

i=l \xi,yi£<8 



Furthermore as the different parts of the p-adic product terms in the definition of C are all 
independent, we can replace all Xi,Xi and yi,y~i by only two variables x,y. By application of 



equation (67) we obtain 



c 'u,v = II X z (.%v,x 1 y)X[(x-y,z-w) u ] 

z, toS© i=l \x,y£® 



E 



X Y\\ ^ Z(u,v,x,y) X [(x y, z w/^l 

i=l \x,y£® 

2p 

^2 Z(u,v,x,y)X[(x-y,z-w) ul ] 
x, y e<$ 



Note that for all z, g G (3 the equation z — w = g has a unique solution w £ 0. Therefore we 
can simplify 

2p 

C ' u ,v = n /^j z^j z { u : v > x iy)% : [( x -y,g)u\ 

gee> x,ye® 



And with the definition of Z(u,v, x,y) as given in equation (66) we get 



gets 



y ] n q ^u-x+y,u-x+y^v— x+y,v— x+y^~ [{ x V: 9)u\ 
x,ye<8 



2p 



6.3 The Affinity Condition (AF) 
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Now, for all x, a G © there is an unique solution y to the equation x — y = a. Thus 

2p 



CL 



n 



4pq- 



g£® 
•,pq+2p+l 

gee 



^ ] ^u-x,u-x^v-x,v-x ' ^ [( x i 9)ui\ 
^1 ^u-x,u-x-^v-x,v-x ' ^ [( X l9)l 



x-e© 



2p 



Defining (7 = n -( 4p9+2p+1 ) • C we have EVAL pin (C) = EVAL pin (C"). We will show that C 
satisfies the statement of the lemma. By condition (R3) we have A q = (3 q + (& q for some 
P q G £5 and (5 g a subgroup of (5. In the following we will make the additional assumption 

that u — /3 q ,v — fi q G Recall that the vertex weights D^\ x U _ X D^ X v _ x are non-zero iff 
!)-i,m-i£A,j. This thus is satisfied iff x G and we may rewrite 



E 
E 

gee 

E 

9 ee> 



2p 



'y ^ D U _ X U _ X D V _ X V _ X • [(x,g) L 
;e©, 

^ [^(u - p q - x) - p q (u - fi q - x) + (g, -x) u ] 
;e© 9 

^2 ^ [P<i( v ~ Pq + X ) ~ Pg( u ~ P g + x) + (g, x) u ] 



2p 



x£®„ 



2p 



where the second equality follows from the group theoretic definition of as given in (R4) 
and from the skew-bilinearity: {x,g) u = — (g,x) w = {g,— x) u . The last equality follows from 
the fact that we are summing over all x G <& q which renders the inversion — x of x irrelevant. n 



We are now in a position to prove the Non-Bipartite Affinity Lemma 6.6 



Proof (of Lemma 6.6). Fix some q G Z w such that q > 0. For every p G N Lemma 6.7 just 
proved shows that the Nymphaea reduction witnesses 

EVAL pin (C) < EVAL pin (F,D) 

for a non-negative real valued nxn matrix C = C(p, q) which satisfies that for all u, v G P q +<3 q 
we have 

2p 



9 e& 



^2 % [Pq( v ~ Pq + X)~ p q {u - P q + x) + (g, X) u ] 



xe©„ 



We need some small simplification and some further preparation before we can give the proof. 
Assume in the following that u, v G P q + <3 q . Consider the term (g, x) u in the above equation, 
where g G (5 and x G <5 q . Define a set K = {g G (3 \ {g, x) u = for all x G & q } which clearly 
is a subgroup of <5 and let &' = &/K be the factor group of & modulo K. Let / : & — > &' be 
the canonical homomorphism defined by g i— > g + K, which thus has kernel ker / = K. Then 
the following holds 
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Claim 1. For all g' G <&' and all a, b G f~ (g') 

we have 



(a, x) u = (b, x) u for all x G & q 



And |<8'| = \<S q \. 



Proof. For the first part, we have a + K = b + K. Thus b = a + k for some k G K which yields 
(b,x) u = (a,x) u + (k,x) u = (a,x) u . Let us now show that \&'\ = |<5 9 |. Note first that the 
columns H*- g = 3£[{*,g)u\ for all g G <£> q form an independent set of columns. Thus there 
is a set I Q & such that the submatrix Hj^ q is non-singular. Note that by the first part of 
the claim, all distinct a,b £ I are contained in distinct fibers of /. That is, there are distinct 
g',g" G <5' such that a G and b G f~ l {g")- This proves = |/| = \<S> q \. H 

By this claim, we may, for each (7 G & , fix an element g G f^ 1 {g) (which is always possible 
as / is surjective) and assume particularly that = 0. We obtain 

2p 



c u ,v — \k\ 



E % [pq{ v ~ Pq + X) ~ p q {u - P q + x) + (g, x)J\ 

xe& q 



Define now C = C(p,q) = \K\^ ■ C. We have EVAL pin (C") = EVAL pin (C) and for u,v G 



Cud — 



^ [pq(v -P q + x)- p q {u - (3 q + x) + (g, x) u ] 



2p 



(68) 



so we may found our further considerations on C' which makes our reasoning a bit more 
convenient. 

Claim 2. For all g G & we have the following: 

(1) If g = then £ [(5, x) w ] = |<5,|. 



(2) If 5 / then 



0. 



In particular, the mappings {x g ■ $ q — > U w | Xg '■= [(g, , g G &} form a basis of the 
\<8 q \ dimensional vector space L 2 {& q ) of functions / : <5 q -»• C. 

Proof. The statement (1) is clear from the definition of For statement (2) recall that <5 q 
has a decomposition into cyclic groups £1 © . . . © (£ z and each x G <3 9 has a representation 
x = J2i=i ^i^i with each /ij being a generator of £j and each Aj G Z ord ( h .y Thus 



a:S© 9 •Vl£ Z ord(hi) ^z£ Z ord(h z ) 



E 



(5>E Ai/i ^° 

i=l 

e n^^^ 



AiSZ ord(fcl) A z eZ ord(hz) i— 1 



= n e ^k^u 

j=l Ai6Z ord(h .) 



A; 



6.3 The Affinity Condition (AF) 
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By g 7^ we know that there is an x £ & q such that (g, x) w ^ 0. Hence by x = Ya=i ^i^i 
there is at least one hi such that (g,hi) u ^ 0. Let N be a multiple of the order of hi. By the 
homomorphism property of the (*, *) u operator we have 



That is X 



% l(g, hiU N = JT [N(g, hi},,} = JT [(g, Nhi) u ] = % [(g, 0) u ] = 1 
is an ord(/ij)-th root of unity an therefore we have 

E %-[(9,hiU Xi = 0- 



/) hi) & 



This finishes the first part of the claim. For the second part note that we have, for g ^ h £ &' , 

= E %~\(9-kx) 



(Xg, Xh) = E XffN • Xftfa) = E ^ ^' ~ x ) 



By the first part of the claim, this value is zero iff g — h £ K which implies g = h by the 
definition of the elements g. Therefore all distinct XgiXh are pairwise orthogonal. The fact 
that they form a basis of L 2 (<5 q ) now follows from |©'| = \<£> q \ as given by Claim[lj _| 

Claim 3. If EVAL pm (ff, D) is not #P-hard then the following is true. For all a £ & q there is 
exactly one g G <5' such that 



\<&a 



^ X [p q {a + x) - p q {x) + (g, x) w ] 



and for all other g' £ &' we have 



^ X [p q (a + x)- p q (x) + (g, x) u 
Proof. As EVAL pin (#,D) not being #P-hard implies the same for EVAL pin (C") we can, for 



all choices of p £ N, derive the following by Lemma 4.23 



For all p £N, the matrix C = C (p, q) does not contain a block of rank at least > 2. 

This is particularly true for all submatrices of C induced by indices from the coset f3 q + <5 q . 
Note that equation (68), Claim [2] and 6 = 0, entail that for all u £ (3 q + & q 



C' 



U,U / j 



E 



2p 



E x [< 6 - 



2p 



\®„\ 2P . 



Consider a 2 x 2 submatrix of C induced by u = (3 q + and some v = f3 q + a £ & q , 



Id* |2p nj 

I I 

u u,i \ VJ q\ 
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As this submatrix is either not part of a single block, or it is of rank at most 1 we see that 

(69) 



For all p £ N either C' UjV = or (C' u J 2 
By the choice of u, v we have 



l<3 9 | 4p - 



^2 & [Pq( a + p q {x) + (g, X)u] 



x£®„ 



and note that the inner term 
For all < i < \<5 q \ define 



2p 



is always at most \& q \. 



g£<5' 



^2 % [Pqi a + X) ~ pq(x) + (g, x)J\ 



This enables us to rewrite the expression for C' u v by 

[« e | 

c^ = Y2caAv 2 T- 

And C' uu = C' vv = \<S q \ 2p . Assume first that C' u v = holds for infinitely many p. Then 
Lemma 4.4 implies that c aj o = |<8'| and c a ,i = . . . = c 0) |g | = 0. However this would mean 
that the mapping / : & q — y U defined by 

f(x) := %~ [p g (a + x) - p g {x)] 

is orthogonal to all Xg f° r 9 G & in contradiction to these being a basis of L 2 (& q ) — as we 
have seen in Claim 2l Therefore, there are infinitely many p such that (C' u v ) 2 = \0 q \ 4p . That 
is, there are infinitely many p such that 

\& q \ 2p = Y2 c ^ 2 r 

u=0 



which by Lemma 4.4 implies that c a ,o = — l and c a ,_ 



C a,|<5,[-1 = and c a) | 09 | = 1. 



This finishes the proof. 



Assume from now on that EVAL pm (i7, D) is not ^P-hard. Then the uniqueness statement 
of Claim [3j implies that there is a mapping 7 9 : (3 q — > such that for all a £ <3 q we have 



I0 O 



% [pq{a + x)- p q (x) - {7/ q (a) t x)u\ 



Further, by 



< 



% [pq(a + x) - p q (x) - (7q(a) , x) u ] 
J2\^'[p q (a + x)-p q (x)-(j q (a),x) UJ ]\ = \<5 q 



x£&„ 



6.4 The Proof of the Non-Bipartite Case Lemma 



6.1 
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we see that equality can be achieved only, if there is a constant a a £ f2 such that for all 



x G ®g we have 



p g (a + x) - pg(x) - (7 g (o), x)u = a a . 
With x = and p 9 (0) = we see that a a = p c (a). This finishes the proof. 



6.4 The Proof of the Non-Bipartite Case Lemma 6.1 



Let H be an w-algebraic n x n matrix and D a family of diagonal matrices defining an 



(H— STD)-problem. Lemma 6.3 implies that EVAL pm (//, D) is #P-hard unless the group 
condition (GC) is satisfied. 

Assume therefore that (GC) is satisfied and recall the definitions of the Representability 



6.4 



Conditions (R1)-(R5) on page [12} We know that H and D satisfy (Rl) and (R2). By 
Lemma 



the problem EVXL pm (H, D) is #P-hard unless the conditions (R3) and (R4) are 



satisfied. By Lemma 6.5 we may further assume w.l.o.g. that (R5) is satisfied. Therefore 



the proof follows as Lemma 6.6 proves #P-hardness of EVAL pin (iT, D) unless the Affinity 
Condition (AF) holds. 

7 Polynomial Time Computable Partition Functions 



In this section we will prove Theorem 3.8 The following two lemmas yield the proof. 



Lemma 7.1. Let H be an u-algebraic nxn matrix and D a family of diagonal matrices defin- 
ing an (H— STD) -problem which has a representation as given in conditions (Rl) through 
(R5). Assume further that the Affinity Condition (AF) is satisfied. Then the problem 
~EVA~L pm (H, D) is polynomial time computable. 

Lemma 7.2. Let A be an co-algebraic and D a family of diagonal matrices defining a (B 
H— STD) -problem which has a representation as given in conditions (B Rl) through (B— 
R5). Assume further that the Affinity Condition (B— AF) is satisfied. Then the problem 
EVAL p '"(yl, D) is polynomial time computable. 



As in the previous sections, we will give the proof only of Lemma 7.1 The proof of Lemma 7.2 
is omitted due to its similarity to the former. For completeness, we refer to [Thu09| for a 
proof. 



The proof of Lemma 7.1 relies on a reduction of the given partition function to functions 
Z q (f) closely related to degree 2 polynomials / over the ring Z 9 (q being a prime power). 
In the first part of this section we will introduce these functions and show how they can be 
computed in polynomial time. Afterwards, we will show how to reduce the computational 
problems satisfying the preconditions of the Lemma to these problems. This will prove both 
results. 

7.1 A Polynomial Time Computable Problem 

Let q be a prime power and / G Z q [Xi, . . . ,X n ] a polynomial. Define Q = exp(2-7ri • q^ 1 ) - 
we fix this definition for the rest of this section. Define the function 

z q {f)= E cp--^. 

Xi,...,X n £Zg 
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We say that / is of degree 2, if there are constants c^-, Cj G Zg such that 

n 

f(Xi ,...,X n ) = J2 CijXiXj + a X i + co- (70) 

i<j i=l 

Let EVAL((7) denote the problem of computing Z g (f) given a degree 2 polynomial / G 
Z 9 [Xi, . . . , X n ]. Since it has been shown in |CCL09] that this problem is polynomial time 
computable, it will be central in the proof of polynomial time computability of partition 
functions. 

Theorem 7.3 (Theorem 12.1 in |CCL 09j). Letq be a prime power. The problem EV AL(q) 
is polynomial time computable. 



7.2 Computing Partition Functions — The Non-Bipartite Case 



We will now prove Lemma 7.1 Recall that we are given an u;-algebraic n x n matrix H and 
a family D of diagonal matrices defining an (H— STD)-problem which has a representation 
as given in conditions (Rl) through (R5) and the Affinity Condition (AF) is satisfied. We 
shall show that the problem EVAL pm (i7, D) is polynomial time computable. 

The proof relies on the reduction of EVAL pm (i7, D) to the problems EVAL(g) which we 
have introduced in the preceding section. To make our partition functions at hand more 
amenable for this reduction, we will first analyze the structure of the vertex weights repre- 
sented by the family D a bit more. 



7.2.1 The Structure of the Mappings p c 

Recall the mappings p c as given in condition (R4) and assume further that they satisfy 
condition (R5) and (AF). In this section we will derive an important property of these 
functions. Before we state it, a word of preparation is in order. Recall that the image of 
all mappings p c lies in the Abelian group Q defined in condition (Rl). We will work with 
binomial coefficients which we define by Q) := x ( x ~^ for all A € Z. Observe in particular that 
(2) is defined also for negative values and for each j£fl the expression L) • g is well-defined. 

Lemma 7.4. For every c G Z w the mapping p c satisfies the following. Let \\,...,X Z G Z 
and h\, . . . ,h z G <3 C then 

\i=l / i=l ^ ' j<i 

Furthermore, for all g G & c the following holds. If ord(g) is odd, then the order of p c {g) 
divides ord(g). If ord(g) is even, then the order of p c (g) divides 2 • ord(g). 

Recall the Affinity Condition (AF): For all c G Z w there is a 7 C : <8 C — > © such that 



Pc{y + x) - p c {x) - p c (y) = {i/c(y),x} u for all x, y G <5 C . 



(71) 



7.2 Computing Partition Functions - 



The Non-Bipartite Case 
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Lemma 7.5. The following holds for all c G Z w and y, y', x £ & c . 

(ic(y + y'),x) u = {ic(y),x) u + {ic{y'),x) u and {j c (x),y)oj = {<y c (y),x) lJj . 



Proof. By commutativity of the addition x + y equation ( 71 ) yields 

(lc{y),x) w = p c (y + x)- p c {x) - p c (y) = (7 C (re), y) u (72) 
By the bilinearity of (*,*) and equation ([72]) we have, for all y,y' G ® c , (7 c (y + y'),x)u = 



{lc{x),y + y')u, = (7c0c),y) w + (7c ( x),y')oj. A second application of ( |72[ ) on the right hand 
side of this yields (7 c (y + y'),x) u = {'y c {y),x) w + (7 c (y'),x) w . □ 



Proof (of Lemma 7.4). We will first prove that, for all h%, . . . , h z £ & c , 

Pci^hA =^2p c (hi) + ^2(-y c (h i ),h j 



(73) 



\i=l / i=l j<i 

The proof is a straightforward induction on z. The case z = 1 is trivially true. For z > 1, 
recall that by equation (71) 

Pc{y + x) = p c (y) + p c (x) + (7 c (x),y) a , for all x,y £ <S C . 

This yields 



'2-1 



'2-1 



Pc ( ^2 hi J + jO c (^) + (j c (h z ), ^2 hi 



vi=l 



vi=l 



vi=l 



2-1 



2-1 



^Pc(^i) + ^2(lc(hi),hj) w + { p c (/i z ) + ^2(-Yc(h z ),h 



i=l 



]<i 



Here, we used the induction hypothesis for p c ^X)i=i hi ) to obtains the second equality. This 
proves (73). We claim that the following is true. 

For all A € Z and ft € <5 C , we have /O c (A/i) = \p c (h) + ( " ) (7 c (/i), (74) 



Before we prove this, let us see how it helps proving the Lemma. By equation (73) we have 

/ z \ 

y^Ajfoj = ^p c (Ai/ij) + ^AiAj(7 c (/ij),/ij) w . 



vi=l 



i=l 



And substituting p c {\hi), for all i G [z], with the expression given by equation (74) yields 
Pc ^2 Xihi 



\i=l / i=l 

just as in the statement of the lemma 



\Pc{hi) + ( 9 * ) (lc{hi), hi),, + ^ A i Aj(7 c (/ii), /ij) 
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The proof of equation (74). For A = this is trivial and for A > recall equation (73). 
We have 



Pc (A/i) 



\i=l / 8=1 ^ ' 



(75) 



This proves equation (74) for A > 0. It remains to prove the case A < 0. Before we do this 
we need to prove the second statement of the Lemma, namely 

Claim 1. Let g G <5 C . If ord(<7) is odd, then the order of p c (g) divides the order of g. If ord(g) 
is even, then the order of p c (g) divides 2 • ord(g). 



Proof. Let N > be the order of g. That is, N ■ g = and equation (|75|) implies 




p c (0) = Pc (N ■ g) = Np c {g) + ( lc (g), g) u . 

If TV is odd then it divides (^) and therefore ( 2 )<? = 0. Hence (^f)(7e(s)> <?)w = ilc(g), (^)d)ui = 
(lc(g),0)u) = and the above equation thus implies = Np c (g). If N is even then the same 
argument using 2N yields the proof. H 



Let us now prove equation (|74j) for A < 0. Define N as the smallest multiple of 2 • ord(/i) such 
that N + X > 0, then Xh = (N + A)/i and p c (A/t) = p c ((iV + X)h). We may thus apply the 
part of equation (74) which has been proved already and obtain 

'N + A 

'A 



Pc ((N + \)h) = (N + X)p c (h) + 



{lc{h),h) L 



*Pc(h)+[ n )Mh),h) 



The second equality follows from Claim [T] 



From Lemma 7.4 we can derive a rather technical result which will become necessary in the 



proof of Lemma 7.1 



Lemma 7.6. Let g G (5 C be an element of order q = 2 k . The order of the element 2 • p c {g) + 
{lc{g)-,g)ui in %w is at most q/2. 



Proof. By Lemma 7.4 we have 

= Pc(0) = p c (q ■ g) = qpc{g) + 



(ic(g),g) 



(76) 



and the order of p c (g) divides 2q. Note that the order of {^ c {g)-,g)ui divides q by q("f c (g), g)u = 

(ic(g),q- g)uj = (7c(s),0) w = 0. 

Assume first that the order of {'y c (g),g)u equals q. Then the element 2 • p c (g) + {'Jc(g),g}u 
has order at most q/2 only if the order of 2p c {g) is q, as well. Assume therefore, for contra- 
diction, that the order of 2p c (g) divides q/2. Equation (76) implies 

V 



o 



(ic(g),g) 



As q — 1 is odd, q does not divide (|) — a contradiction. 

Assume now that the order of (j c (g)j g)w is strictly smaller than q, that is q/2- {j c (g),g) 
0. As q/2 divides (t) equation (76) implies that = qp c (g)- This finishes the proof. 



7.2 Computing Partition Functions - 
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7.2.2 The Final Reduction 



We will now prove Lemma 7.1 We start with some preliminaries. Consider Zg-polynomial / 
of degree 2 as in equation (70). With each variable Xi we associate an order ord(Aj) which 
divides q. A degree 2 monomial CijXiXj for i ^ j is called consistent if the values -ord(Xj) 
and C{j ■ ord(Xj) are divisible by q. Similarly, a monomial CjAj is consistent if q divides 
Cj • ovd(Xi). Finally, a square monomial caXf is consistent if 2 • cu • ord(Xj) is divisible by q. 

Lemma 7.7. Let q be a prime power. The following problem can be reduced to EVAL(g) 
in polynomial time. The input consists of a degree 2 polynomial f £ X q [Xi, . . . ,X n ] of the 
form given by equation (70). Each variable Xi has an order ord(Xi) which divides q and all 
non-constant monomials are consistent. The task is to compute the value 

E ••• E c^'-' Xn) - 

Proof. To perform the reduction of the problem at hand to KVAL(q) we will need to show 
how to "lift" the order of each variable X, to q. Assume w.l.o.g. that the order x% := ord(Aj) 
of each variable is greater than 1. Fix some variable X v . We will have a look at the situation 
after the order of X v has been lifted. For all % £ [n] define % = q ■ x^ 1 . We have 

E E Cq {Xi ''"' x "''"' Xn) = E EE ( fiXi '-' x " +j ' x "'-' Xn) . (77) 

Let £(Xi, ... , X u , ... , X n ) = Y,i <u c iyXi + ^2 u<i c ui Xi and fix the degree 2 polynomial h 
which captures the part of / not involving X u . We have 

f(X\, . . . , X n ) = c vu Xl + X v ■ £(Xi, . . . , X u , . . . , X n ) + h(X\, . . . , X u , . . . , X n ). 

Further 

f(X\, . . . ,X„ + j ■ x u , . . . ,X n ) = c vv (X v + j ■ x u ) 2 + (X v + j ■ x v ) ■ i + h 

= c vv Xl + 2c vv x v ■ jX u + c vv x 2 v ■ j 2 + X v i + jx u £ + h 

By the consistency of the monomials qXj and that of c^XiXj for all i / j we see that all 
coefficients of divisible by q. This is also true for 2c uu x u . To see this for c vv x u note 

first that if q is a power of 2 then c uu x 2 is divisible by 2c uu x u as 2 divides x u . Otherwise the 
condition that 2c vv x v be divisible by q is equivalent to c vv x v being divisible by q. Therefore, 
for all j e {0, . . . , q u - 1} we have f(Xi, . . . , X v + j ■ x v , . . . , X n ) = c vv Xl + X v £ + h (mod q). 
Using equation ((77]), we see that 



E E cf Vi v - n,- E E cr--^. 

Mi^y. Xi& n x„&z q Xiez^ x n & Xn 

As this holds independently for all v E [n], it follows that 

z q (f(x 1 ,...,x n ))=(fi q )j y: ■■■ e ^■■■■ Xn) . 

\i=i / X!& x , x n ez x „ 

This finishes the proof. 
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7.2.3 The proof of Lemma [7TT 



Proof (of Lemma\7l\). Let now EVAL pin (i?, D) be an (H-STD) problem which meets all 
conditions (R1)-(R4) and (AF). Let G = (V, E) be a digraph an 4> a pinning. With 

V = V \ def(0) we have 



(78) 



0Co-:y-S>© UV&E 



We shall prove that Zh,t>(<I>,G) can be computed in polynomial time. To do this, we will 
present a polynomial time reduction to the problems EVAL(g). Most of this proof will be 
devoted to rewriting the expression of this partition function into a form which resembles the 
expression Z q (f) used in EVAL(g). 



Recall the representation of H and D as given in conditions (R1)-(R5) on page 12 Note 



that a configuration a which maps some vertex v E V' to cr(v) Anj^j does not contribute 
to this partition function. Assume therefore that cr(v) £ A^g^u is satisfied for all v £ V. 
Slightly abusing notation we will denote Ajq^)] by A^ = /3 V + & v . We introduce, for each 
v £ V a variable a v £ & v . For all v E def(0) we assume that the value of a v is fixed to 
a v = 4>(v). Application of the representation (R1)-(R5) thus yields 



EE*" 



v<=V 

^ {Pu + o- u , p v + a v ) u + P(d( v )i (& v ) 

.uv&E v£V' 



Let us rewrite the expression E^ ul)g£; (/3 u + a u , f3 v + cr v ) w a little bit. Denote by /i uv for all 
u, v E V the multiplicity of the edge uv in E. Then, using the skew-bilinearity of the 
operator, 



uvGE 



U,VdV 



For each v £ V define the value % = E^ ug y(/x ut , — fJ, vu )Pu- Then 

uv£E u,v£V v£V 

Let /?„ denote the mapping /0[a(„)] and define a value 



c= 



u,v£V vGde{(4>) 



The value c is constant for all different values the o~ v corresponding to v £ V' may assume. 
We have 
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which, by a straightforward computation yields 



We will proceed by further rewriting the terms of the expression 



(79) 



V>uv{<Tu,o v )u + E ( 7 ?i" cr f)w + Pv{°v)- (80) 

For some technical reason which will become clear later, we now decompose 2u into prime 
powers. Let wi • • • u ( = 2w be a decomposition of 2u into powers of pairwise distinct primes. 
Denote by p v the prime corresponding to oj v and assume w.l.o.g. that uj\ is a power of 2 - 
if u) is not divisible by 2 then we let oj\ = 2. 



Recall once more the representation given by (R1)-(R5) on page 12 We have (8 = 
%q 1 © • • ■ ®^q z for some prime powers It will be convenient to write this decomposition of 
(5 in terms of the above prime decomposition of 2uj. Define sets P u = {i G [z] \ p v divides q{\ 
for each v G [t\. By condition (R2), the order of each g G (5 divides to implying that the 
Pi , . . . , Pt form a partition (with possibly empty parts) of [z] . 

Therefore, we rewrite the decomposition of (5 by © = ©„ =1 ©j g p„ 2^ • Further, with each 
<S C being a subgroup of we see that they also have decompositions (5 C = ©* =1 ©j £ p ^q c < 
such that w.l.o.g. g Cj j divides % for all c£Z u and all i £ [z]. Let g c ^ be a generator of Z 9c . 
for all c G Z w , G [i] and % G P,. For notational convenience, for every v G V we denote by 
g^i the element for all i G [2]. 

Define a vector X = (-X7«,i))' U ey,ie[z]- Every <j v G <5„ has a representation in terms of a 
tuple (X( u> i), . . . , X( U)je )) of variables each satisfying Xr Vji \ G Z Qv i such that 

z 
i=l 

Define matrices T, A and a vector A, with indices in V x [z] such that, for all (u,i), (v,j) G 
V x [z], 

a / {lv(9v,i),9v,j)u, , ifu = veV'andj<i , . 

A WX^) " \ , otherwise. {82) 

E(u,i)(?;j) = A 4 ™ ' {du,ii 9v,j) uj (83) 

A _ / Pv(9v,i) + (Vv,9v,i)u, , if f G V',i G [z] , . 

™ \ , otherwise. 1 J 

Further, define vectors B(X) and H such that for all (v,i) G V x [z], 



/A 



2 "> l) ) , if u G V",* G [z] 



B1X >- - r ' ' 11 (85 » 

(86) 



(lv(gv,i),9v,i)w , if u G V, i G [z] 
"^'^ 1 , otherwise. 
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We can now rewrite the terms of equation (80). Applying equation (81) and the definition of 
r in (83), we find that 



u,»£V u,v£V i,j=l 

Further we have, for each v G V' 



xrx' 



i=l 
z 



vi=l 



J2 x (v,i) (pv{gv,i) + (vv,gv,i)u) + (^^J ^ v (Sv,i),9v,' 



j<i 



where the second equality follows from Lemma 7.4 By the definition of A, A, B(X) and 5 
this can be rewritten into 



^2(rh,<Tv) u + Pv(<rv) = X • A T + XAX T + B(X 



Define /(X) = XrX T + X • A T + XAX T + B(X) • E T . By equation ((79} we therefore have 

Z H) u(<f>,G) = c-J2^[f(^)} (87) 
x 

where the sum is over all vectors X = (-XV^, i))( v £\ 6 y x rj such that Xr v q G ^q vt for all (v,i) G 
V' x [z\. For all (v,i) G def(</>) x [z] the values X/ V n are fixed so as to satisfy a v = (f>(v). 
Note that this is always possible by the definition of a v in equation (81). We call any vector 
X satisfying these conditions a valid assignment. 

We will simplify the description of / further. Define X = {I„ | v G [t]} with /„ = | 
v G V, i G Pj,} — That is, the sets 7^, form a partition of V x [z]. Recall that Xj and Ajj 
denote subvectors and submatrices of X and A induced by the indices in I and J. Then 

/(X) = Y, X 7 r 7J Xj + X.AjjXj + Y, X/ • Aj + B(X)/ • 3? 

/,Je2 /ex 

We will claim that the "cross-terms" in this expression are zero, i.e. 

/(X) = Y Xjrj/Xf + XjAjjXf + Xj • Aj + B(X)/ • Ej. 

To see this, we shall prove that, for all I, J £ I with / / J we have Tjj = and Ajj = 0. 
We show this for Ajj, the result follows analogously for T. By definition (cf. equation (82)), 
for v G V, i G I and j G J we have A^v^-) = {^ v {9v,i)-,9v,j)ui- By the properties of 7„ (cf. 
Lemma 7.5) and the bilinearity of {*,*)aj we have N^y^gyj), g v ,j)ui = v,i)i 9v,j)u) — 

(lv(gv,i),Ng Vj j} i j. Therefore the order of A( U)i )( v j) divides both, ord e (g V}i ) and ordg^j). 
By our choice of I and J these orders are co-prime and the claim follows. Define 

/„(X/) = X/r 77 Xj + X/Az/Xj + X 7 • Aj + B(X)/ • 3j for every i/ and J = I„. 
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Thus 

t 

/(X)= J>(X/„; 



Claim 1. The order of all entries of T, A and S divides u and the order of all entries of A 
divides 2ui. 



Proof. For A, T and H this follows directly from their definitions in equations (82), (83) and 



(86). Every non-zero entry of A satisfies A(„ j) = p v (g v ,i) + (Vvi9v,i)u} and the claim follows 
from Lemma |7.4| -I 



Claim 2. Let v G [t] and I = I u , then all non-zero entries of the matrices Tu, Ajj and the 
vectors A/ and 5/ are divisible by oj^ for all fi ^ v. 

Proof. For (u,i), (v,j) G I u we have ^(u,i),(vj) = u. uv ■ (g u ,i, 9v,j)u- The order of this element 
(in Z2 W ) in is some divisor of co u and therefore the element itself is divisible by 2u ■ (tjj v ) . 
The proofs of the other statements follow analogously. H 

Claim [l] implies that /(X) G for all valid assignments of X. Define ^2w( x ) = exp( 2 ^ x ). 
We see that 

ZhM, g) = c J2 W • K x : /( x ) = a ( mod 2uj )}\ 

By Claim [2] we see further that /(X) = a (mod 2ui) is equivalent to a system of simultaneous 
congruences given by f u (X.i v ) = a (mod uj v ) for all v G [t]. We have 

t 

Zh,t>{4>,G) = c- [a] |{X/„ : f u (X Iv ) = a (mod ui u )}\ 

By the Chinese Remainder Theorem each a can be written uniquely as a = Ylu=i uf' ' a v f° r 
values a u G [0, u) v — 1]. Recall that by Q q we denote the value ( q = exp(27ri • Then 



r 2oo 



a h 



|{X/„ : /i,(X/J = a v (mod uj u )}\ 



zhm,g) = c - e ••• e 

"lG^^ OL t &L ut u=l 
t 

= ° n e c • k x /» : = <*» ( m ° d ^)>i 

The partition function Zh,t>(<P,G) thus simplifies to 

and its computation reduces to the computation of t many independent values of the form 



E^: (X/J - (89) 
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We shall show now, that each of these problems satisfies the preconditions of Lemma 7.7 
This then implies that each problem is reducible to EVAL(a; i/ ) and therefore polynomial time 



computable by Theorem 7.3 



Fix some v E [£]. To show that the problem (89) satisfies the preconditions of Lemma 7.7 



we have to show that the monomials in /„ are consistent. In particular, the square monomials 
have to be consistent. Recall that this means that they are of the form b • X? such that 



2 • b ■ ord(Xr V)i -\) = 2 • b ■ q v ^ is divisible by q. By equation (88) and I = I u , we have 



/„(Xj) = X/IYrXj + X/Aj/Xj + X 7 • Aj + B(X)/ • Ej . 

Observe that the condition a v = 4>(v) for all v G def(0) turns some variables X( V j) into 
constants. Therefore, some degree 2 monomials turn into linear ones and others of these and 
some linear ones turn into constants. Therefore, to show that the consistency preconditions 



of Lemma 7.7 are satisfied, it suffices to do this for the original polynomial /„. For the 
terms X/F/jXj + X/A//XJ the consistency of the corresponding monomial follows from 
their definition. It remains to consider the monomials arising from X/ • + B(X)j • Ej . 
Recall that 



X 7 • Aj + B(X)/ • Ej = ^2 X (v,'i) • (Pv(9v,i) + (Vv, 9v,i)uj) + ( X( 2' 4) ) {lv{9v,i),9v,i 



(90) 



(v,i)el 



Case A. v / 1. We have u u a prime power of some odd prime p u . Note that 2 has a 
multiplicative inverse in Z Wj/ and therefore the term 



X 



2' l) )(lv(9v,i),9v,' 



X (v,i) \ X (v,i) - 1) 



\'yv\9v,i)i 9v,i> j: 



inequation (90) produces consistent monomials. Further, the term X( v ^y{p v {g v ^) + {r] v , g Vt i)uj) 
is consistent as well. To see this, note that Xr v ,i) ' {Vv, 9v,i)u is consistent by definition of Xr V i) 
since ord(Xr Vti \) = ord(g Vt i). By Lemma 7.4 the term Xr v ^ ■ p v {g Vi i) is consistent as the order 
of g v i is odd. 



Case B. v = 1. By equation (|90|) we have 
X 7 • Aj + B(X)/ • Ej = 



^ X( Vji ) ■ (pv(9v,i) + (r] v ,gv,i)uj ~ 2 1 ■ {lfv(9v,i),9v,i 
(v,i)ei 



" 2 • 2 1 • {'7v(gv,i),9v 



Note that the expression 2 1 • (^v(gv,i), gv,i)u is well-defined, as (jv{gv,i), gv,i)uj has, by defi- 
nition of u)\, order at most uj\/2. This proves that for the above square terms we have 

2 • ord(X (V)i) ) • 2~ l ■ (j v (g v ,i),g v ,i) u = ord(X (V)i) ) ■ (nf v (g v ,i),g v ,i) u 

is divisible by u\. Further, X( v a ■ (r] v ,g Vt i) w is consistent. And it remains to show this for 
X (v,i) ■ (Pv(gv,i) ~ 2 _1 • ( r Yv(9v,i),9v,i)u) — which follows from Lemma 



7.6 
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